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ABSTl^CT 

In  the  previous  report^^\  ^e  fundamental  equation  that 
describes  limit  cycles  in  nonlinear  sampled-data  systems  has 
been  derived.  In  that  case,  the  equivalence  of  limit  cycles  with 
finite  pulsed  systems  havii^  a  periodically  varying  sampling-rate 
I  is  observed,  and  the  methods  of  analysis  applied  to  the  latter  are 

extended  to  obtain  these  limit  cycles  with  the  aid  of  final  value 
theorem. 

This  fundamental  equation  is  modified  and  is  simplified  to 
some  extent  under  certain  assumptions  as  it  can  be  applied  to 
systems  both  with  and  without  integrators.  The  limitation  on 
the  longest  period  of  saturated  and  unsaturated  oscillation  is  in¬ 
vestigated  and  the  critical  gain  for  their  existence  is  derived, 

,  starting  from  the  modified  fundamental  equation.  Also,  the  stability 

of  limit  cycles  as  well  as  the  equilibrium  point  is  considered,  based 
on  Neace's  method  and  its  modification. 

t 

Through  this  study,  various  kinds  of  non-linearities,  namely, 
pulse -width  modulation,  relay  saturating  amplifier  with  linear  sone 
and  quantised  level  amplifier  are  discussed,  and  examples  are 
presented  for  each  of  these  cases.  Self-excited  oscillation i  are 
mainly  examined,  as  well  as  the  possible  existence  and  stability  of 
limit  cycles,  however,  the  method  can  be  extended  to  forced 
oscillations. 

Finally,  experimental  works  are  performed  in  order  to  verify 
the  theoretical  results  by  means  of  digital  computors. 
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CHAPTER  1 


Introduction 

Nonlinear  sampled-data  feedback  control  systems  have  been 
extensively  investigated  in  recent  years  (Fig.  1).  One  of  the  most 
important  topics  of  these  studies  has  been  the  stability  problem  of 
such  nonlinear  sampled-data  systems.  Lyapunov's  second  method 
has  been  introduced  as  a  powerful  tool  for  the  analysis  of  the  stability 
problems.  ^  At  the  same  time,  the  classical  describing-function 
method  has  been  used  in  analogy  with  the  case  of  nonlinear  continuous 
systems.  ^ 

Among  the  various  types  of  nonlinearities  the  pulse-width  modu- 

g 

lation  has  been  studied  by  such  people  as  R.  F.  Nease  ,  R.  E. 

Andeen?'  E.  Polak,  “  T.  T.  Kadota,  ^  W.  L.  Nelson,  Shao  Da 

Chuan,  I.  V.  Pyshkin,  and  E.  I.  Jury  and  T.  Nishimura^' 

9  10 

For  example,  R.  E.  Andeen  *  presented  an  approximate  method  for 

the  analysis  by  replacing  the  pulse-width  modulator  by  the  equivalent 

pulse-amplitude  modulator.  Hence,  his  method  is  restricted  by  the 

small-signal  condition  because  of  that  approximation.  On  the  other 

hand,  Shao  Da  Chuan^^  gave  an  exact  analysis  of  the  limit  cycle  in  the 

PWM  system  using  the  canonical-form  expression.  However,  his 

analysis  is  limited  to  the  limit  cycle  of  two  sampling  periods.  I.  V. 

14 

Pyshkin  presented  an  extensive  method  to  prove  the  existence  of 
limit  cycles  in  PWM  systems,  but  still  he  had  to  resort  to  certain 
approximation  using  the  low-pass  characteristics  of  the  linear  plant 
and  the  stability  problem  was  not  substantially  attacked. 

T.  T.  Kadota^  applied  the  second  method  of  Lyapunov  for  the 
stability  discussion  and  succeeded  in  deriving  a  sufficient  condition  for 
asymptotic  stability  in  the  large  of  PWM  systems,  although  his  condition 
was  rather  conservative  compared  with  the  actual  stability  condition 
obtained  experimentally. 

E.  I.  Jury  and  T.  Nishimura^'  derived  a  method  to  find  the 
exact  behavior  of  limit  cycles  in  PWM  systems  extending  the  theory 
of  the  periodically-varying  finite-; pulse-width  systems.  Also,  they 
discussed  certain  features  of  the  stability  of  such  systems  in  connection 
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FIG.  I  NON-LINEAR  SAMPLED  DATA  CONTROL  SYSTEM 
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with  the  stability  of  limit  cycles  when  they  are  in  the  simplest  form, 
namely,  of  the  two  sampling  periods. 

In  this  report,  the  method  introduced  in  the  above  references 
is  extended  and  simplified  in  its  form.  Hence,  this  report  is  in 
direct  continuation  with  the  previous  one  and  further  detail  of 
that  report  will  be  summarized  in  the  last  part  of  this  introduction. 

The  main  subject  of  this  report  is  the  stability  problem  of 
PWM  sampled-data  systems.  By  the  Lyapunov -function  method  one 
attempts  to  find  the  Lyapunov  function  which  is  positive  definite  and 
is  decreasing  for  every  sampling  instant.  In  order  to  find  such  a 
general  function  one  has  to  expect  the  worst  case  might  happen, 
without  paying  attention  to  the  various  features  of  disturbances  within 
the  system.  This  fact  leads  to  the  result  that  the  sufficient  condition 
derived  by  this  method  is  often  excessively  conservative,  as  is  seen 
in  the  case  of  Ref.  6. 

In  contrast  to  such  a  comprehensive,  macroscopic  method,  we 
adopt  the  microscopic  method  in  which  we  start  by  sepaiating  the 
possible  modes  of  oscillation  into  the  adequate  categories,  then 
analyse  the  stability  of  oscillations  of  each  category  from  its  simplest 
mode  and  forward  the  analysis  until  we  cover  all  the  possible  modes 
of  oscillations. 

This  laborious  work  may  appear  almost  prohibitive  in  its 
beginning.  However,  by  finding  the  certain  regularities  which  govern 
such  oscillations,  we  can  achieve  the  precise  analysis  of  the  stability 
problem.  This  is  the  basic  attitude  of  this  report  towards  the 
specified  problem. 

In  Chapter  II  the  simplified  form  which  yields  the  feature  of 
limit  cycles  will  be  derived  and  this  equation  will  be  extended  to  include 
the  system  which  contains  pure  integrators. 

In  Chapter  III  we  will  derive  the  method  to  find  the  longest  period 
of  limit  cycle  of  relay  mode  oscillation.  Then  we  will  discuss  the 
stability  boundary  for  the  relay  mode  oscillation.  The  same  technique 
will  be  extended  to  unsaturated  oscillation  as  well  as  to  other  nonlinear 
systems. 

In  Chapter  IV  the  stability  of  the  equilibrium  point  will  be  in¬ 
vestigated  and  final  conclusions  will  be  derived  on  the  stability  of  the 
PWM  sampled-data  systems. 


Review  of  the  Previous  Report^ 

We  present  a  short  review  about  the  content  of  the  previous  report 
entitled  "On  the  Periodic  Modes  of  Oscillations  in  Pulse-Width  Modulat¬ 
ed  Systems"^  and  reproduce  those  equations  which  will  be  used  in  this 
report.  The  scheme  of  the  PWM  system  is  shown  in  Fig.  2.  The  PWM 
controller  (lead  type)  has  such  characteristics  that  its  output  is  a  unit 
pulse  (positive  or  negative),  the  sign  of  the  pulse  is  identical  with  the 
sign  of  the  control  error  at  the  sampling  instant  (e^(0))  and  the  pulse 
width  is  proportional  to  the  magnitude  of  e^(O). 

Since  the  sampling  period  T  is  fixed,  the  maximum  pulse  width 
is  T  and  the  saturation  will  occur  beyond  this  point.  Then  the  output 
of  the  PWM  controller  during  the  (n-fl)th  sampling  period  will  be  given 
by  the  product  of  the  sign  function  ‘Y(n)  and  the  unit  pulse  function  ^(t). 
The  origin  of  the  time  axis  t  is  placed  at  the  nth  sampling  instant. 


where 

en(0) 

"Te;TUTT 

and  this  takes  only  -l-l  or  -1. 


(l.l) 


(1.2) 


Uh(t)  =  1 


for  0  <  t  <  h 
—  —  n 


=  0 


elsewhere 


(1.  3) 


Also  the  pulse  width  is  given  as  follows; 


=  «|e„(0)| 


for  0  ^  a  I  ®^(0)  <  T 


=  T 


for  a  I  e^(0)  ]  >  T 


(1.4) 


where  a  is  a  gain  of  the  PWM  controller.  The  incremental  response 

AC^(s)  is  defined  as  the  product  of  the  transfer  function  KG(8)  and  the 

Laplace  transform  form  of  the  pulsed  output  e',  (f)  from  the  PWM 

nn 

controller 

AC^(s)  =  KG(8)E;^(b)  (1.  5) 
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nT)i-e(T) 


PWM  6|,(T)  linear!  c(r) 
CONTROLLER  PLANT 


FIG.  2o  PULSE  WIDTH  MODULATED  SYSTEM 


{n+2)T  (n+3)T 


(n*2)T  (n*3)T 


nT  (nti)T  |»J  L — J  (n+4)T  (n+5)T  i 


e(T)  ,  e*(T)  ;  INPUT  a  SAMPLED  INPUT 

TO  CONTROLLER 

eJ,(T)  ;  OUTPUT  OF  CONTROLLER 
FIG.  2b  INPUT  a  OUTPUT  OF  PWM  CONTROLLER 


where 


-h  s 


(1.6) 


I'lu.'*'!  =  V(n) 


Then  the  incremental  response  of  each  limit  cycle  is  given  by  the 
summation  of  AC^(8)  during  one  limit  cycle,  multiplied  by  the  relative 
delay  factor 

M-1 

Z^C^(s)  =  ^  AC^(s)e"*®'^  (1.7) 

1=0 

The  z-transform  of  this  becomes 

M-1 

AC*(a)=  ^  **^ACj(*)  (1.8) 

1=0 

Applying  the  skip-sampling^^  operation  to  AC*(z)  yields  the  incremental 
response  at  the  instants  0,  T^,  ....  where  T^  =  MT,  the  period 

of  the  limit  cycle; 

ACg^(Z)=  ZjACJ(z)]  (1.9) 

Then  it  is  shown  that  the  response  at  the  beginning  of  the  limit  cycle  is 

obtained  by  adding  AC'*'(Z)  for  all  the  periods  of  limit  cycles  and  by 

*  17 

applying  the  final  value  theorem  to  the  summation.  In  this  case,  it 
is  assumed  that  KG(8)  has  no  integrator  in  order  to  insure  the  convergence 
of  limiting  process. 

Thus  the  final  result  is  derived  as  follows; 


c 


Os 


lim  AC*  (Z) 
Z^l 


=  lim 
Z-»l 


.-I 


> 


[KG(s)E'j^(8)] 


(1. 10) 


This  equation  can  be  applied  for  any  type  of  nonlinearities.  In  the  case 
of  the  PWM  system,  E|^(a)  in  the  above  equation  is  replaced  by  Eq.  (1.  6). 
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Ex.  (1.  10)  can  be  extended  in  general  form  to  give  the  response  at  the 
jth  sampling  instant  of  the  limit  cycle. 


c.  =  lim  (Z) 

Z-^1 


s  lim  Z 
Z-^1  • 


M-1 

^  z"^[KG(e  E'^j^(a)] 


since 


♦Z) 


Z^[z^  AC*  (z)] 


(1.11) 


(1.12) 


IVhen  KG(s)  has  an  integrator,  the  problem  is  solved  with  the  aid  of 
difference  equations. 

Using  these  equations,  the  limiting  cycle  of  M  x  2  and  M  x  4  are 
analyzed  in  the  examples  of  the  previous  report^  and  the  existence  of 
limit  cycles  of  the  PWM  mode  and  the  relay  mode  as  well  as  the  stable 
region  is  indicated.  Also  the  critical  gains  for  each  region  are  specified 
for  the  limit  cycle  of  M  x  2. 

In  the  appendix,  the  stability  of  the  limit  cycle  is  discussed,  which 
is  another  important  problem  of  the  limit  cycle  as  much  as  its  existence. 

g 

The  discussion  is  based  on  the  important  theorems  given  by  Nease. 

The  definition  of  the  stable  limit  cycle  is  given  as  follows. 

DEFINITION:  A  limit  cycle  is  said  to  be  absolutely  stable  if  any 
small  perturbation  about  the  limit  cycle  approaches  0  as  n  oo. 

Then  two  theorems  are  referred  to  which  present  the  method  for 
testing  the  stability  of  limit  cycles. 

THEOREM  3.  Assume  that  the  nonlinear  difference  equations 

have  a  periodic  solutions  S^]  of  period  M,  and  that  the  functions  F(X  ) 

are  single  valued  and  possess  continuous  first  partial  derivatives.  The 

first  approximation  of  the  difference  equation  for  small  perturbations 

about  this  periodic  solution  S  1  is 

n'' 
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I 


i 


and  the  solution  of  this  equation  determines  the  stability  of  the  periodic 
solutions  S^]  if  [  nonsingular  at  all  of  the  solution  points.  In  the 

above  equation  is  the  perturbation  about  S^]  and  the  components  of 
[  A  1  are 


a..  =  3F.(X  )/8X., 

ij.n  i'  n'  j 


X  =  S 
n  n 


(1.15) 


THEOREM  4.  The  system  of  £q.  (1. 12)  is  stable  if  all  the  eigen¬ 
values  of  the  matrix 

lie  inside  the  unit  circle.  Then  all  the  solutions  tend  to  0  as  n  becomes 
large. 

Theorem  4  of  Ref.  8  is  partially  stated  in  the  above,  considering 
only  the  necessary  part  pertinent  to  this  discussion.  Combining  these 
two  theorems,  the  condition  for  the  stability  of  limit  cycles  will  be 
reduced  to  the  following  statement  that  "all  the  eigenvalues  of  the 
matrix  [A^]  which  consist  of  the  first  partial  derivatives  of 
F(X^)  lie  inside  the  unit  circle  at  all  the  periodic  solution  points 
S  ].  "  With  this  condition  satisfied,  the  small  perturbation  will  tend 
to  0  as  n 00,  hence  the  limit  cycle  is  said  to  be  stable  according 
to  the  definition  of  previously  defined  stability  of  limit  cycles. 
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CHAPTER  11 


Modification  of  rundamoatal  Equation  for  Limit  Cycles 


The  general  form  o£  the  aquations  which  gWe  the  raspoasa  at  tha  Jth 
sampling  instant  of  the  limit  cycle  in  tha  nonlinear  samplad-data  system  has 
been  derived  in  Eq.  (1. 11  U  is  assdmed  that  KO(s)  has  no  poles  at  the  origin. 

Certain  modification  of  Eq.  (i.  lli)  is  attempted  in  the  following  part. 
Since  the  content  of  the  skip-sampling  operator  [  ]  is  a  function  of  s,  i.  e. , 
a  sampled  function  with  the  period  T,  the  operator  Z^  sam|des  tha  sampled 
function  of  T  with  the  period  T^  ■  MT.  Hence  this  samptiag  process  is  redun¬ 
dant  and  the  first  sampling  process  with  period  T  may  he  removed.  For  that 
purpose,  the  content  of  the  Z  operator,  i.  e. ,  ficf  (s)  shall  be  modified 
in  such  a  way  that  ACj^(a)  contains  the  incremental  response  of  one  period 
of  limit  cycle  preceding  the  Jth  sampling  instant.  Hence,  the  new  incremental 
response  AC*^(a)  is  described  as  follows. 

r' 

l«0 


i  •  0,  1 . M-1. 


It  is  understood  that  the  system  has  already  been  on  the  limit  cycle,  hence 
has  the  periodical  feature  with  the  period  T^. 

wli«l.l>M  (.mi 

Substituting  AC^^(b)  for  AC*^(b)  in  Eq.  ( 1- 11  )  yields  the  desired  response  at  the 
ith  sampling  instant. 


c 


is  • 


llm  Z  lAC*(a)J 
Z-ai  • 


r  M-1 


lim 

Z-H 


I 

l«0 


»*'  y  iKO( 


1  .  0,  1,  . . 


M-1 


Inverting  the  content  of  the  operator  to  the  function  of  s,  and  replacing  the 
skip  sampling  operator  Z  by  the  ordinary  L  operator  with  the  period 

■  W  %dmT* 

the  above  equation  becomes  as  Eq.  (ird  ),  given  below,  where  Z  ■  •**** 


! 


ii  multiplied  to  the  content  of  the  brocket  while  Z“^  le  multiplied  to  the  Outolde 
of  the  bracket. 


c.  «  lim 
Z  — 


.1 


(•“"K 


eT 


lim 

2-e 

s  -I 


1 

0 


z 

HT 


-1 


I 


e^P^KO(p) 


.4" 


1.0 


dp 


i  -  0.  I, 


M-l 


U4  ) 

MeT 

The  multiplication  of  e  to  the  content  of  the  bracket  ie  done  to  oaeure  the 
convergence  of  the  integrand  to  aero  along  the  Infialte  aomlcircle  on  the  right« 
half  plane,  and  thie  doee  not  cauae  any  change  in  the  final  reeult  under  the 
condition  that  Z  — ♦  1. 

The  limiting  proceee  e  0  may  be  performed  before  the  integratioa  if 
the  path  of  integration  T  along  the  jw  aaie  ia  taken  eulficiently  cleee  to  the 
imaginary  axia  of  the  p« plane. 


Hence,  Eq.  (2.  4  )  ie  reduced  to 
M-1 


ie 


1 

ZwT 


I 


.“"■flKilrt 


2 


I  •  • 


tmO 

UpT 


dp 


Thie  is  the  modified  form  of  Eq.  (2*  1  ). 

Further  aimplification  is  poesible  when  the  oscillation  ie  aymmetrical  and 
monotonic.  The  symmetrical  oecttlation  is  such  that  the  same  ahapo  of  oeciUa* 
tioa  is  repeated  for  every  half  period  of  the  limit  cycle  with  an  opposite  sign. 
The  word  monotenic  in  this  case  implies  that  the  input  SigU)  lo  1^*  nonlinear 
component  is  positive  for  the  first  hall  period  of  the  limit  cycle  and  ia  negative 
for  the  other  half  period.  These  conditions  are  esproesod  mathematically 
as  follows: 

Letting  M  ■  2p  (p  «  number  of  samples  in  half  period  of  limit  cycle) 
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{ l.f>) 


ft 


for  I  ■  0,  1, 


and 

y(I  )  m  ^  \  for  I  ■  0,  1,  ... ,  •! 

(2.  7  ) 

Y<l)  •  -1  for  I  •  K.  . 

Tkt  akovt  attumptloa  it  Juatifiod  ia  mott  of  tht  catot  ia  practlct  bocaatt  the 
liatar  plaat  usually  potttttst  tht  low*patt  eharaettriotict,  which  rsadtr 
tht  otcillatioat  smooth  aad  moaotoaic. 

Alto  it  it  tmphatiatd  that  tht  attumptioa  of  tymmttry  it  plactd  oa  tht 
iaput  of  tht  liatar  plaat,  at  ia  Bq.  aot  oa  its  output.  Htaet,  tht  atym- 

metrical  otcillatioa  with  rttptct  to  tht  plaat  output,  at  obttrved  ia  tht  ttcoad- 
ordar  relay  tytttm,  (Fig.  9'oi  reference  1)  is  still  the  symmetrical  oscillation 
■  with  respect  to  the  plant  input  and  satisfits  the  condition  of£q  (  2.  6) 

Siact 


Eq.  (  5 )  it  rowrittaa  at  follows 


(  2.  8) 


c 


it 


1 


I 


77^ - 


••HT 


J* 


777** 


(2-9  ) 


That,  half  of  the  tarot  of  the  dtaomiaator  of  the  iattgraad,  which  iacludtt  the 
taro  at  the  origia,  art  caactlltd  by  tht  atrot  of  (  1  <■  ia  the  aumarator 

by  iatrodttciag  the  ceaditioa  of  ayauattry  aad  moaoteatcity  of 
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oscillations.  ^  This  ellmiaation  of  the  pole  at  the  origin  yields  a  great  contri> 
bution  for  the  analysis  of  systems  with  integrators  which  shall  be  developed 
in  the  following  part. 

Equation  (2. 10)  will  be  rewritten  for  the  case  when  the  function  of  the 
nonlinear  component  is  the  pulse-width  modulation.  Then, 


-  Yd) 


(2.11  ) 


where  y(I)  i*  specified  by  £q.  ( 1.  2  ,  and 


h 


I 


Substituting 


c 


is  " 


a|e^(0)| 

«  I  r|(0)  -  c^(0)| 


Eq.  (  2. 11 )  into  Eq.  (  2. 10)  yields 


(  2. 12) 


dp  (2. 13) 


So  far,  it  is  assumed  that  KG(s)  has  no  integrator.  The  problem  of  KG(s) 


The  conditions  of  symmetry  and  monotonicity  are  the  sufficient  coaditions  for 
cancelling  the  pole  at  the  origin  ol  the  integrand  of  Eq.  (2. '06),  but  are  not  the 
necessary  coaditions.  Actually,  the  condition  of  monotonicity  is  not  required  at 
all  in  order  to  cancel  the  pole  at  the  origin.  Only  the  condition  that  every  *11.(0) 
in  one  limit  cycle  have  its  pair  of  the  same  magnitude,  but  of  opposite  sign,* 
in  the  same  limit  cycle,  is  required  for  that  purpose.  In  special  cases,  this  condi¬ 
tion  is  further  reduced  to  that  the  summation  of  el.  (0)  is  equal  to  sero,  as  in  the 
case  of  the  quaaitaed  level  amplifier. 

For  example,  the  four-period  limit  cycle  of  1,  k,  -k,**l,  |  k|  <1,  which  is 
shown  in  the  example  of  Ref.  19 for  the  saturating  amplifier  system,  satisfied  the 
above  condition. 

However,  this  does  not  assure  the  existence  of  such  limit  cycle  as  shown 
in  the  same  reference  in  which  the  existence  of  the  above  limit  cycle  is  denied. 
Moreover,  the  stability  of  such  asymmetrical,  non-monotonic  limit  cycle  is 
quite  dubious,  as  is  noticed  in  the  experimental  observation.  Therefore,  the 
introduction  of  the  assumption  of  symmetry  and  monotonicity  is  practical  when 
the  linear  plant  has  the  low-pass  characteristics,  as  in  most  casss.  Also,  we 
emphasise  that  this  assumption  does  not  introduce  any  approximation  for  the 
analysis  at  all. 
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with  integrators  will  be  discussed  in  ths  following  part. 

An  example  of  second-order  system  with  a  single  integrator  is  sotred  in 
the  previous  section  with  the  aid  of  differsaco  oquntlono.  The  reason  for  the 
difference  equations  being  used  in  addition  to  our  fundamental  equation  le  that 
the  pole  at  Z  «  1  which  originates  from  the  single  integrator  of  the  ptaat  aaahos  the 
application  of  the  final  value  theorem  impossible,  aince  when  Z  approacheo  unity 
the  term  AC*(Z)  would  diverge  because  of  the  pole  at  Z  ■  1.  This  fact  can  alee 
be  observed  when  we  take  a  look  at  the  modified  equation  (S«4  ). 

When  the  integral  is  evaluated  by  the  residue  method  for  all  die  poles  of 
KG(p),  the  term  1  -  causee  the  integrand  to  diverge  when  the  residue 

at  p  ■  0  is  evaluated,  and  when  s  —a  0  and  ^is  dees  net  give  any  finite  value  for 
the  residue. 

However,  if  the  condition  of  symmetry  and  monotonicity  is  introduced  to  such 
a  system,  this  troublesome  pole  at  the  origin  is  eliminated  as  s  — ^  0  and  the 
evaluation  of  the  residue  at  the  origin  becomes  possible.  This  is  seen  in  Sq. 

( e«10)  ^  which  the  pole  at  the  origin  of  the  integrand  has  already  been  removed. 

.  Therefore  KG<P)  nay  contain  not  only  a  single  integrator,  but  also  a 
double  or  triple  or  any  higher  order  of  integrator  at  the  origiiL  And  it  is 
proved  that  the  simplified  equation  ( Z.  10)  is  valid  for  any  shape  of  tte  glMt 
KG(s)  if  it  satisfies  the  physically  realisable  condition  under  the  assumption  sf ’ 
symmetrical  and  monotonic  oscillation.  Also.  Eq.  (2. 10)  is  applicable  for  any 
type  of  nealinearities  if  the  output  af  nonlinear  element  can  be  specified 

as  a  function  of  the  input  and  output  of  the  over-all  system. 

We  have  derived  three  fundamental  equations  for  limit  cycles  in  this 
chapter,  as  shown  in  Eqs.  (l'U>).  (  Z.  5).  and  (  Z.  9  ).  Wo  will  explain  briefly 
the  advantage  and  dieadvantage  of  using  each  of  these  equations. 

The  first  equation  of  Eq.  ( 1.  U  )  has  ths  disadvantage  of  having  two  kinds 
of  s -transformation,  namely,  the  a-transformation  wifii  respect  to  T  and  the 
skip-sampling  operation  with  respect  to  T^.  However,  in  actual  calculation,  we 
can  perform  the  skip  sampling  operation  by  picking  up  the  necessary  terms  out 
of  the  expansion  in  powers  of  s*^,  without  carrying  over  the  integration  of  sampl- 

^WOthse.  fioreed  to  nse‘J^V'4^UK/whfp  system 

gijmnlai  n  fanctidn  o£ «.  not  o£  s. 
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The  eccond  equation  of  Eq.  ( 2.  b  )  may  be  used  in  moot  caooo  except 
when  the  oyetem  hao  digital  procesolng  unite.  However,  we  muet  be  careful 
in  carrying  over  the  integration  becauee  the  integrand  coataine  the  delay  factor e 
in  the  form  of  e**^^  and  more  labor  ie  required  for  the  computation. 

The  third  equation  of  £q.  (2.  9  t)  is  ueeful  when  the  oecillation  ie  eymmetric 
and  monotonic,  and  ie  the  moot  convenient  form  among  theee  fundamental 
equatione.  When  theee  aeeumptione  are  violated,  ae  obaerved  in  the  examplee 
of  the  quanlteed  level  amplifier  of  Section  2.  2,  (3),  we  have  to  uee  either 
Eq.  (1.  11)  or  Eq.  (2.  5  ). 

When  the  plant  hae  an  integrator,  Eq.  ( 2.  5  i)  or  Eq.  (  2.  9  )  can  bo  uaed, 
although  they  give  only  the  oacillatione  of  aero  D.  C.  component.  ^  In  caae  the 
oecillation  of  non-eero  D.  C.  component  ie  deeirad  or  the  valuea  of  derivativee  of 
reaponaea  are  required,  the  difference  equatione  are  very  helpful  aa  ehown  in 

section  V  of  refcrcnct?  1. 

Example  2.  1 

The  equatione  which  yield  the  eolution  for  thia  eymmetrical  PWM  mode 
oecillation  are  illuatrated  in  thie  example. 

The  plant  ie  choaen  aa  second  order  and  the  period  of  the  limit  cycle  ie 
M  r  2p  sampling  periods. 

In  Eq.  (2. 13),  l«t  i  ■  0  and  take  ytifi)  ■  y(| )  ■  'fl  for  all  I,  referring  to 
Eq.  {,2.  7  ).  Hence  Eq.  (2.13)  becomes  as  follows. 


O.  C.  component  implies  the  average  value  of  reeponees  at  the  sampling 
instants  during  one  period  of  limit  cycle. 
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'Oi 


I  I 

nr  I 


p(p+B)  C 


dp 


1-0 


K  f  “o  .-^■’^(.“’o-i.  .  ■H  . 

tg/.-HW,'  *  TT-  ^  - 


h  ,  .-bT,  ’*(1-1  ,, 

^  bd  ♦ 


I 


(  2. 15) 


Multiplying  -  o  to  both  aides  of  the  above  equation  and  knowing  hQ  «  -a 
the  following  equation  is  derived: 


h 


0 


bh. 


(2.16) 


For  1  f  1,  It  Is  observed  from  Eq.  (  ^  )  that 

y(l+l)-+l  for  I  -0,1 . 1a-2 

«  •!  for  I  •  p>l 

and 


Hence,  h^(s  obtained  in  a  similar  manner 


(2.17) 


b,  .  -j-i 


[L  [■’  J 


•[ 


(  2. 18) 


In  general. 


for  f  -  0,  1 . p>l*l 

for  I  ■  p«l,  p'lfl,  . . .  ,  p-1 


-  >1 
a  •! 


(2. 19) 


and 


*V-l+l  ■  N-l 

Than,  uaing  thaaa  raaulta 


h 


i 


Ml  + 


•] 


Ml  ♦  a^^) 


U.20» 


i  ■  0,  1,  . . . ,  |i*l 


Tha  aolutiona  ol  tha  aquations  u4ieh  ara  glvan  by  Bf.  (2.  20)  far  I  ■  0,  1,  ...  fk*l 
yiald  tha  axact  faatura  of  tha  PWM  moda  oocUlatioa  that  wUl  bo  aaatalaad 
within  tha  cloaad-loop  PWM  aamplad-data  ayatoaa.  Allhao^  dio  aamarlcal 
aolutiona  of  thaaa  traaacandaatal  aquatiooa  ara  not  oaatty  obtaiaod  u4aa  p  >  2, 
cartain  approximatioaa  for  tha  axpoaaaiial  tarma  ara  poaaibta,  aa  (Uacaaaad  ha 
tha  prariotts  aaetian.  Alao  tha  programming  on  tha  high*apaad  digital  aaaapatar 
will  anabla  ua  to  aolva  thaaa  aquatiena  without  appronimatian  triMn  tim  oyatam 
conatant,  tta  aampliag  pariod,  and  tha  oyatam  gaiaa  ara  apaalflod. 
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Example  2.  2  Quantised  Level  Amplifier 


111  tkls  SHaaipts,  ths  awdi— r  gala  amplifisr  has  a  charactsristic 
of  qoaatissd  Isvol  (2  Isvsls)  as  shoiwa  ia  Fig.  (3). 

This  charactsristic  is  rsprsssatsd  as  follows. 


.;(t)  .  i 

■  0.  s 

•  0 

•  -0.  5 
«  *1 


if  d  <  s JO) 

m 

if  d/2  <  sJO)  <  d 
■ 

if  •d/2.<  sJO)  <  d/2  (  2.  21 ) 

if  -d  <  sJO)  <  -d/2 

tt 

a  sjo)  <  -d 

B 


whore  ad  ■  1 

The  liaear  plaat  is  agaia  secoad  order  givea  hy  E%.  (2. 14). 

We  will  deiBoastrate  the  existence  of  asyauaetrical  oscillalioa  with  the 
period  of  three  sasapliag  periods. 

For  such  aa  asyounetrical  oscillatioa,  Sq.  ( 2.  5  )  mast  be  used. 

Letting  M  •  3,  and  it  is  assumed  that  the  output  of  the  honliaear  compoaeat  has 
a  sequence  of  0,  -1,  heace  qfO)  ■  1,  q(2)  ■  -1.  Substituting  these  values 
into  Eq.  (2.  5  ),  together  with  Eq.  (2.21  .),  yields 


os 


'is 


K  2T 

¥ 


1 . 


K  ft  I  -  1 

*  P  'm'.. J 


] 


(2.  22  ) 


(  2.  23  ) 


(2.  24>) 


The  stahiltty  of  this  limit  cycle  is  tested  by  the  method  which  is 
described  ia  Chapter  Xi  as  well  as  by  the  digital  computer  oxperimeat. 

It  is  fouad  that  ^s  limit  cycle  is  unstable  and  will  move  to  the  relay  mode 
oscillatioa  of  M  ■  4. 
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FIG.  3  CHARACTERISTIC  OF  QUANTIZED  LEVEL  AMPLIFIER 
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When  T  «  2,  b  ■  1,  K  ■  1  and  d  ■  1/3  (a* 3),  those  values  are  calculated 
as  follows. 

c  .  .  0. 

os 


®2s 


0.  565 


The  output  of  the  quantised  level  amplifier  corresponding  to  c^^  is  0 
for  c^^,  -1  torc2^,  satisfy  the  original  assumption  of  the  sequence. 

This  asymmetrical  limit  cycle  is  shown  in  Fig.  ( 4 ). 

Another  example  of  an  asymmetrical  limit  cycle  will  be  demonstrated 
for  the  same  system  as  the  previous  example,  except  that  d  ■  1/2  (aaZ) 
in  this  case.  The  sequence  of  the  output  from  the  nonlinear  component  is 
assumed  as  >1.  0,  -1/2,  0.  1/2,  0,  -1,  having  seven  samples  in  one  period 
of  limit  cycle. 

Again  £q.  ( 2.  5  )  is  used  for  the  analysis  of  such  asymmetrical 
oscillation.  For  example,  the  response  at  the  first  sampling  instant  of  the 
limit  cycle  is  given  by 


os 


K 


5 

7 


«  -0.  571 


-bT. 


,.^>T  ,  ..4I.T  . 


■] 


2b(l  .  e 


(2.21,) 


Similarly,  the  responses  at  the  other  sampling  instants  of  the  limit  cycle 
are  obtained  as  follows. 


is  • 

-0.178 

c,  «  0.  470 

2s 

Cj^  ■  -0.  011 

4S  * 

-0.  371 

'5.  ■  0-'^’ 

'6s  •  ® 

Applying  these  values  of  responses  to  E^.  (2.,  21  .)  (letting  e JO)  ■  -  cJO)  ), 
we  find  that  the  assumed  sequence  will  be  reproduced  from  the  nonlinear 
gain  amplifier  and  such  a  limit  cycle  will  be  maintained.  The  feature  of 
this  limit  cycle  is  shown  in  Fig.  ( 5<).  These  two  limit  cycles  are  proved  to 
be  stable  limit  cycles  by  the  digital  computer  experiment. 
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1 


C| 


®n  '  1 1  0,  - 1 


FIG.  4  ASYMMETRICAL  LIMIT  CYCLE  OF  QUANTIZED 

LEVEL  AMPLIFIER  SYSTEM 
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1 


c 


C| 


i  +  1 , 0, -0.5 , 0,+0.5, 0, -I 


FIG.  5  ASYMMETRICAL  LIMIT  CYCLE  OF  QUANTIZED 

LEVEL  AMPLIFIER  SYSTEM 
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CHAPTER  in 


ANALYSIS  or  SATURATED  AND  UNSATURATED  OSCILLATIONS 
AND  THEIR  STABILITY  BOUNDARIES 

In  th«  previous  chaptort  wo  have  dorivod  a  fuadamoatal  oquatioa  ta 
trace  the  exact  behavior  of  limit  cycles  which  ars  sustaiaod  withia  the 
nonlinear  sampled-data  systems.  We  have  aloe  domonstratsd  that  such 
oscillations  can  be  eliminated  by  reducing  the  gain  in  PWld  systems  when 
they  are  in  the  simplest  mode.  However,  the  problem  of  finding  gain 
boundary  to  eliminate  all  the  possible  limit  cycles  requires  an  enormous  amount 
of  work. 

But  we  will  proceed  to  solve  the  stability  problem  In  this  chapter  by 
the  localised  approach  to  saturated  and  unsaturated  osciUntlens  in  PWlf 
systems. 

3. 1  Limitation  on  the  Period  of  Limit  Cycles  of  Relay  Idede  Oscillations 

It  has  been  observed  that  there  exists  a  certain  limitatien  on  the 
longest  period  for  the  limit  cycle  which  Is  sustained  within  the  autonsmens 
relay  sampled-data  system.  Inawa  and  Weaver^^  discussed  this  problem 
on  the  second- order  system  and  derived  an  equation  u^ch  gives  the  maxi- 
mxim  half  period  of  limit  cycle  as  a  function  of  the  sampling  period,  intro¬ 
ducing  a  fictitious  delay  to  the  sampler.  Also  Pyshkln^^  treated  the  same 
problem  using  the  dencribiag-functien  method. 

We  will  show  that  this  problem  can  be  solved  for  any  order  of  die 
linear  plant  without  any  approximation.  Under  the  assumpCien  of  syssmetric 
and  monotonic  oscillation  at  the  output  of  the  nonlinear  component,  udilch 
is  an  ideal  relay  in  this  case,  Eq.  (2. 13)  is  valid  for  the  relay  mode  eocillatien 
if  all  h|^^  are  replaced  by  the  sampling  period  T.  Hence 

•i.  •  •  I  ^ 

r  ‘  foO 


22 


It  it  to  W  rocolUd  tiMt  tk«  output  ot  tko  Mollaoor  olooooat  lo 

pooitlTo  for  tho  first  half  poriod  aad  aogoitvo  for  tko  otkor  kalf  porlod. 

This  imptiot  tkot  tho  roopoooo  C|^(0)  !•  Mgatlro  for  tho  first  half  poriM 
and  positiro  for  tho  othor  half  if  tho  aoaliaoar  eompoaoal  such  as  rolaf^ 
PWhl,  or  saturatiag  amptifior  has  ao  hystorosis  ia  it  (ao  iaput  is  assupsod). 
Thoroforo, 

y(l)  ■  >1  for  I  ■  0,  1 . p«l 

■  -1  for  I  ■  Pt  P^l . M-1 


and 


'Is 


>  0 


for  I  ■  0,  1 . p*l 

for  I  ■  p,  P'flt  . . .  •  M-1 


<3.  3> 


Thoa  tho  maaimum  aumbor  p^^  of  oaaiptiag  poriodo  ia  tho  half  poriod  of 
limit  cyclo  is  oasUy  ohtaiaod  hy  Urostigatiaf  tho  polarity  of  ^  , 

which  is  tho  rospoaso  at  tho  last  sampliag  iastaat  ceaiaiaod  ia  tho  first  half 
poriod  of  tho  limit  cyclo,  aad  is  givoa  hy  Ef .  (3. 1),  lottiag  i  ■  p*l  ia  it^ 


Vi.»  *  '1^  J 

r 


ff5»“  “ 


yflfp»l)o"*'^  dp 
(3.|4) 


Siaco  is  aogativo  hy  tho  original  assumption,  tho  polarity  of  c^_^^  ^ 
is  tostod  starting  from  p  >  2,  i.  e. ,  from  c^^  .  ^  may  ho 

aegatiro  up  to  a  cortain  p  ■  p',  but  may  hocome  positiro  for  all  p  >  p'-fl. 

In  that  caoo  this  critical  p*  is  takon  as  p^^^^  siaco  g  must  ho 
aogativo  by  tho  original  assumption,  aad  whoaoror  this  assumption  is 
▼iolatod  wo  may  coacludo  that  such  limit  cyclo  with  tho  half  poriod  pT 
eaaaot  ho  sustaiaod  ia  that  systom. 

This  coaclttsioa  will  ho  illuotratod  for  first«ordor  aad  socoad*opdor 
plants  ia  Eaamplos  3. 1  aad  ).  2. 


Eaamplo  3.1 


Whoa  tho  plaat  traaofor  fuactioa  is  first  ordor,  giroa  hy 


(3.5) 


KG<t)  .  ^ 


The  respoaee  ^  at  the  margin  of  the  flret  half  wave  ia  aimilarly 
obtained.  Letting  i  •  p-l  in  Eq.  (2. 19  ) 


and 


■  fl 

■  *1 


Cqi  Cj,  ....  c  .  <  0 


for  I  ■  0 

for  f  «  1,  2,  .  . . ,  !»•! 


(3.  M 


Hence  Eq.  (3.  4)  becomea 


.  K  .  . . ,  +  ..-(•‘-UbT,  „  g, 

I  ^  e 


It  ia  clearly  obaerved  that  c^^^  ^  ia  poaitlve  for  all  p  >  2,  violating  the 
original  aaaumption  of  Eq.  (3.  7).  Therefore  it  ia  concluded  that  the  firat* 
order  relay  aampled-data  ayatem  can  maintain  only  the  relay  mode  oacillation 
of  two  aampling  perioda  (pal)  and  that  the  longer  period  oacillation  cannot 
be  auatained  ia  that  ayatem.  The  above  atatement  ia  true  for  any  plant  gain 
K,  plant  ^time-conatant  b,  and  aampling  period  T. 

Example  3.  2 

Next  the  caae  of  the  aecond« order  ayatem  with  an  integrator  will 
be  diacuaaed.  The  plant  tranafer* function  ia  given  by  Eq.  (2. 19  ). 

Then,  letting  all  ^  in  Eq.  (2. 13)  of  Example  2.  1  yielda  a 

aet  of  aolutioaa  c^^,  c^^ . ^  for  the  reaponaea  at  all  aampling 

inatanta. 
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'i.  ■  -  ^ 


w 

For  ^  ■  1  *ad  i  ■  «  0 


(.4.  L^’  W«*>>  J 

■  A-.  I-  ’ '  w  *  1 


(3.9) 


®0« 


K  Pr  1  - 


bTO  ♦ 

I  ♦  •* 


(3. 10) 

It  can  be  oatily  proved  that  the  numerator  i&eide  the  bracket  ie  poeitive  for 
all  bT  — ^  0.  Hence  It  hae  been  ehown  that  c^^^  le  negative  for  all  T. 

Then  for  pel  and  i  e  p-1  ■  1 


c 


le 


.  .-«,2 

TTTsf- 


<  0 


(3. 11) 


Therefore,  c^^  ie  negative  for  any  eampling  period  T. 
Next  for  p  e  3  and  i  e  p>l  e  2  in  Eq.  (3.  9) 


®2e 


K 

¥ 


JLl* 

Ml  t-  c 


4-  e 


-bT 


(3.U) 


Simple  trial  will  ahow  that  ie  negative  for  email  T,  but  when  T  becomee 

larger  C2^  will  become  poeitive,  violating  the  original  aeeumptien.  Thua 

we  can  find  a  certain  critical  T^  bu  ept^ing  Eq.  (3. 12)  for  T  when  C2^  in 

equated  to  aero.  Then  the  oecillation  ofpa  SorMeO  can  eidet  for 

T  <  T^  but  it  cannot  be  auetained  for  T  >  T  . 
c  c 

A  nimilar  phenomenon  ie  obeerved  for  larger  p.  Hare  we  may  derive 
an  equation  which  yielda  the  relatien  between  the  iwbebrfw**!"*  p  and  the 
eampling  period  T.  ItOttiag  i  ■  p-1  in  Eq.  (3. 9), 
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K 

1? 


(3.13) 


K  f  T  .  ..bT 


^.(M-l-l)bT  ,  ^-bT, 

b(l  + 


( 


bT ,  bT 

•  2)  ~  6 


,-t.bT 


777^ 


•f  • 


bT 


} 


By  the  original  aeeumption,  ^  <  0.  Hence  ie  obtained  by 

finding  the  maximum  ^  which  eatiefiee  the  following  inequality 


2)-  e 


bT 


1  -  .-•‘W 

777^ 


■f  e 


bT 


-  1  <  0 


(3.14) 


Thue  we  are  led  to  the  following  interesting  conclusion,  that  the  limit  cycle 
of  two  sampling  periods  and  four  sampling  periods  can  exist  in  the  relay 
sampled«data  system  for  any  sampling  period  T,  but  that  for  p  >  3  the 
maximum  number  of  sampling  periods  which  can  be  contained  in  one  limit 
cycle  is  restricted  by  certain  conditions  that  are  specified  by  the  system 
constant  b  and  sampling  period  T,  as  shown  in  Eq.  (3. 14).  m  a 

function  of  the  sampling  period  T  is  plotted  on  Fig.  (•). 

The  discussion  and  conclusion  so  far  can  directly  be  applied  to  the 
relay  mode  (saturated)  oscillation  in  PWM  systems.  The  fact  that  there 
exists  a  limit  on  the  period  of  the  limit  cycle  for  the  relay  mode  oscillation 
must  be  remembered  and  this  will  be  very  useful  in  stability  discussion  of 
PWM  systems  that  will  be  given  later. 

3.  2  Gain  Boundary  of  Relay  Mode  Oscillations  in  PWM  Systems 

It  is  observed  in  Section  IV  of  reference  1  that  the  oscillation 
will  be  shifted  from  the  PWM  mode  oscillation  to  the  relay  mode 
oscillation  when  the  over>all  gain  aK  is  increased,  because  of  the 
saturating  characteristics  of  pulse  width  modulator.  The  critical 
gain  between  these  two.  modes  is  denoted  as  aK^  and  is  obtained  for 
the  limit  cycle  of  two  sampling  periods  in  Section  IV  of  reference  1. 

It  is  not  a  difficult  problem  to  obtain  the  critical  gain  for  longer  periods, 
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FIG.6  MAXIMUM  NUMBER  OF  SAMPLING  PERIODS  CONTAINED  IN 

HALF  PERIOD  OF  LIMIT  CYCLE 


■lace  wc  hair*  already  darirad  Iha^aaaral  aquatioa  for  tha  PWM  made 
oaclllatioa  in  tha  praaleua  aaetloa}  Tha  IraMltlea  will  occur  uAaa  all 
pulaai  h|  of  aaeh  aampliag  period  reach  aaturatloo  aad  fill  every  aampliag 
period  of  limit  cycle.  However*  we  have  derived  the  eqtuitieo  of  c^^  ae  a 
function  of  a  K  in  Sq.  (2.13)  and  we  know  for  the  autooemoue  ayatem  that 

“i  • -'i.  . 


when  the  period  of  the  limit  cycle  ia  choaea  aa  J)tT.  Then  the  critical 
gain  a  ie  given  aa  auch  a  gain  that  every  h|^  becemee  aimultaaeoualy 
equal  to  T  on  both  eidea  of  Eq.  (2, 13);  thue  the  aum  of  h|  for  the  half  period 
becomea  equal  to  the  half  period  pT. 

Combining  Eqo.  (2. 13)  and  (3. 15)  and  lettiag  h|~-->  T  yielda 


pT 


lim  y  h.  ■  lim  (  »a  /  c.  \ 

•>.— f  So  \  So  “j 


- 

•  Tvl 


1^1 

i  I 


oK  r 

nr  J 


e^**^0(p)  ^  y<ltU*****^U  - 
1-0 


UO  r  (1  ♦  e****^)p  ^ 


dp 


0(pMl  -  e 


\  UO  tmO  / 


(3. 16) 


dp 


Taking  out  the  termf  under  the  double  eommatlon  iaaide  the  bracket  of  the 
integrand  and  denoting  thie  ae  f(p,  pT),  we  have 

L-l  J»-1 


f(p 


i«o  i«o 


U^l  U»1 

.  ^  >-l-l)pT  2  . 

leO  UO 


(S.17) 
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Actual  Talus*  of  ^l^i)  aro  flvon  by  Eq.  (  2.  W).  Subotltutlag  their  Taluce 
into  the  abore  yield*  the  general  form  of  f(|»,  pT)  ae  follows 

PT)  .  p  .  2  ♦  (p  .  4)***’’  +  (p  -  ♦  . . .  -  (p  - 


(3. 18) 


Denoting  the  gain  which  satisfies  Eq.  (3. 15)  as  aK^(p)  it  is  giTsn  by 


It  must  he  pointed  out  that  the  assumption  of  the  elmultaneous 
saturation  of  the  pulee  width  h^  is  not  practical  except  when  pal,  since 
the  ware  chap*  of  the  response  may  become  somewhat  like  a  sine  wave, 
but  may  not  appear  like  a  square  wave.  Therefore  most  of  o  |  c^^|  becomes 
larger  than  T  when  all  the  pulse  width*  are  saturated,  so  the  left  hand  side 
of  Eq.  (3.  lb)  may  be  larger  than  pT. 

Hence,  aK^(p)  which  is  derived  from  Eq.  (3.  lb)  must  be  Megarded 
as  the  lowest  gain  that  might  allow  the  existence  of  the  specified  saturated 
oscillation;  it  give*  the  sufficient  condition  for  the  non*exist*nce  of  the 
specified  saturated  osclUation.  Also  it  may  be  regarded  ae  the  critical 
gain  when  the  square  wave  assumption  is  introduced  as  explained  later,  then 
it  gives  approximately  the  gain  when  the  avergge  of  a  |  c^^l  is  equal  to  T. 

The  above  equation,  Eq.  (3.19),  shall  be  applied  for  the  case*  of  the  first- 
and  second-order  systems  in  Examples  3.  3  nnd  3.  4. 

Example  3.  3 

In  the  first  place,  when  the  transfer  function  of  the  plant  is  the  first 
order,  given  by  Eq.  (3.  5),  and  whenp  ■  1,  oK^  is  easily  calculated  as 
follows: 

^(P.  pT)  •  -1  forpel  (3.20) 

Hence,  by  Eq.  (3. 19) 
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aK 


1 

liT 


I 


T 

yr 


bT 


(1  +  eP^)p 


T 


i  >  •-■’T 
- :rST 


(3.  21) 


dp 


Example  3.  4 

When  the  plant  it  second  order,  given  by  £q.  (  2. 14),  the  integral 
of  Eq.  {Z  19)  is  first  evaluated  as  follows 


1 


r  (1  -  eP^lfGt,  pT) 


p  (p+b) 


j^hpT  + 


■bT) 


1  4  e 


(3.  22) 


where  %,  >bT)  is  obuined  by  substituting  •bT  for  p  in  Eq.  (3.18).  Substituting 
this  into  Eq.  (3. 19)  yields  a  as  a  function  of  p; 


aK,(|i)  . 


- il>VT»  .  .•“"I  _ 

b^TU  +  e***'*^)  +  2(1  -  f  (p,  -bT) 


(3.  23) 


a  is  actually  calculated  for  p  «  1,  2,  3,  4  as  follows 

2b^T(l  4  e'**"^) 

aK-fl)  ■  '  '  *  uif  -a—..-  y- 

bT(l  4  e  ®^)  -  2(1  -  e'*^^) 


(3.  24) 


bT - BT- 

•  bT(l  4  e*^**^)  -  e''*T(l  -  e'*»T) 


(3.  25) 


_K<31  .  _ |b^T(l  4  e*^^"^) _ 

•  3bT(l  4  e'^**^)  4  2(1  -  e’^Ml  -  -  Se’***^ 


) 


(3.  26) 


■*r4\)T  ^  '  '~»8bT~'  '  3kT'  '  26a) 

•  2bT(l  4  e  '‘*^)  4  2(1  -  e  '*^)(1  -  e  -  Ze'^**^) 
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It  was  found  in  the  previous  section  that  s  limit  ealsts  for  the  longest 
period  of  the  relay  mode  oscillation  in  the  PWM  syotom  as  well  as  in  the 
relay  system. 

When  this  restriction  is  combined  with  the  results  about  the  critical 
gain  for  the  relay  mode  oacillation  we  reach  an  important  conclusion  con¬ 
cerning  the  stability  boundary  of  the  relay  mode  oscillation. 

Since  we  know  how  to  find  u_^.  for  the  relay  mode  oscillation  and 

nimx 

also  the  boundary  gain  corresponding  to  each  p,  wo  can  find  out  the  lowest 
gain  for  any  relay  mode  oscillation  by  inspecting  every  e  as  a  function  of 
p,  where  p  ranges  from  1  to  p^^-  We  may  conclude  that  no  relay  mode 
oscillation  can  exist  below  that  gain,  since  we  have  covered  all  the  possible 
modes  of  relay  oscillation. 

This  lowest  gain  min(a  K^)  is  the  stability  boundary  with  respect 
to  the  relay  mode  oscillation,  and  it  is  clearly  a  stifficient  condition  for 
the  non-existence  of  that  type  of  oscillation. 

This  derivation  of  the  stability  boundary  for  the  relay  mode  oscillation 
shall  be  illustrated  on  the  first-order  and  Second-order  plants  in  the 
following  part. 

In  case  of  the  first-order  plant,  it  has  been  wfiown  that  only  liaait 
cycle  of  two  sampling  periods,  he.,  pal,  can  exist,  and  we  derived  the 
critical  gain  o  for  this  mode  in  Eq.  (3.  21).  Hence  this  aK^  is  the 
stability  boundary  for  the  relay  mode  oscillation.  When  the  gain  is  lower 
than  this  a  the  system  is  stable  as  far  as  the  relay  mode  oscillation 

is  concerned.  However,  it  is  pointed  out  in  Chapter  n  that  this  oK  enactly 

a  * 

coincides  with  the  abeOlute  stability  boundary  derived  by  Kadota  using 
Lyapunov's  second  method. 

Therefore,  the  question  will  arise  whether  wo  may  take  the  gain 
boundary  for  relay  mode  oacillation  as  the  absolute  stability  beuadary 
or  not,  and  if  the  answer  is  affirmative,  how  we  can  justify  it.  This  question 
will  be  answered  in  the  following  part. 

Next,  the  second-order  plant,  will  bq  discuesed.  In  this  esmo,  the 
problem  is  not  so  simple  as  Ih  the  first-order  syetem,  since  we  iMwi  Itat 
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th«  limit  cycla*  of  two  and  four  oampliag  periods  always  escist  for  any  T 
and  the  longer  period  oscillation  can  exist  when  T  becomes  shorter.  However, 
we  can  prove  fortunately  that  oK^(p)  is  the  monotonically  decreasing 
function  of  p,  hence  the  lowest  gain  boundary  can  be  obtained  only  by  calculating 
o  for  which  is  specified  by  Eq.  (3. 14)  for  a  given  T.  The  proof  of 

the  decreasing  characteristic  of  oK^  is  given  in  Appendix  B.  Thus,  the 
steps  to  follow  in  finding  the  stability  boundary  are  first,  find  by 

£q.  (3. 14)  for  a  specified  T,  then  find  oK^  corresponding  to  this  by 

Eq.  (3.  24);  ultimately  this  o  stability  boundary  for  the 

relay  mode  oscillation  for  that  T.  The  curve  of  this  boundary  is  plotted 
in  Fig.  7  . 

3.  3  Gain  Boundary  of  Unsaturated  Oscillations 

In  the  previous  sections,  we  have  derived  the  sufficient  condition  for 
non-existence  of  relay  (saturated)  mode  oscillation.  Extension  . -of  the  same 
technique  to  the  unsaturated  oscillation  is  attempted  in  this  section.  The 
problem  is  not  so  simple  in  this  case  because  the  pulse  width  may  take  any 
value  between  sero  and  T  and  is  not  fixed  at  T  as  in  the  saturated  oscillation. 
The  derivation  of  the  equation  which  gives  the  longest  period  for  the  unsatu¬ 
rated  oscillation  as  well  as  the  equation  which  yields  the  boundary  gain  for 
such  oscillation  may  become  prohibitively  complicated,  since  each  pulse 
width  at  each  sampling  instant  differs  from  the  other  and  they  cannot  be 
easily  calculated. 

However,  by  introducing  the  sine  wave  approximation  or  square  wave 
approximation  we  can  eliminate  such  complexity  to  some  extent  and  can 
follow  almost  the  same  steps  as  in  the  case  of  saturated  oscillations  to  reach 
the  gain  boundary  of  unsaturated  oscillation. 

We  consider  the  square  wave  approximation  in  which  we  assume  that 
the  pulse  widths  are  constant  and  are  equal  to  h,  which  may  be  considered 
as  the  average  value  of  pulse  widths  during  one  period  of  limit  cycle. 

In  that  case,  letting 
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GAIN  BOUNDARY -FOR  SATURATED  OSCILLATION  AND  STABLE  REGION  OF  PWM 

OSCILLATION  OF  M  =  2 


in  Eq-  (2. 13)  yi^id* 


•  TiT 


r  *  ♦  •  imo 


(3.  *7) 


Then,  we  can  find  which  givci  the  longeet  period  admitted  for  the 

uneaturated  oscillation  by  testing  the  polarity  of  c  ,  .  aad  using 

®  OIMi 

The  gain  boundary  is  given  as  follows  by  modifying  Eg.  O* 


oK^(p.h) 


1 

JiT 


J  fil£lS 


(3.2t) 


The  applications  of  these  two  equations  shall  be  illustrated  on  the  first  order 
and  second  order  systems. 

Example  3.  5 

In  case  of  the  first  order  system,  the  equation  to  give  e  ,  becomes, 


P-l.s 


K 


(3.  29) 


1  e 


Observing  Eq.  (3.  29),  the  same  conclusion  is  derived  as  in  the  saturated 
oscillation  that  the  longest  period  of  the  unsaturated  oscillation  is  2T,  and 
the  oscillation  of  p  >  2  cannot  be  sustained  in  the  first  order  system  with 
the  square  wave  approximation. 

Hence  we  find,  for  uneaturated  oscillation,  that 


u  <* 
'^max 


I 


The  corresponding  gain  can  be  easily  calculated  from  Eg.  (3.  2i)  as  foUaws 


oK,(h)  * 


oK. 


(3.  90) 
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This  coincldet  with  the  eritle«t  galB  of  th«  PWI4  oacUlatioa  aa  darivad  la  dia 

Section  IV  of  reference  1 


Eatampie  6 


The  aaeond  order  ayatem  la  dlaciaaaad.  Under  the  aquara  wave 
aaaumptloa,  the  equation  for  c  .  la  eaaily  derived  from  Bq.  (3.  27)  aa 
fellowa 


c 


P-1.» 


.•*-1  I  - 


•f  a 


hh 


(3.  31) 


Then  the  aaaxlmum  number  of  aamplea  which  can  he  contained  Ini  one  half 
of  the  period  of  limit  cycle  of  unaaturated  oaclUatlon  la  obtained  hy  finding 
the  maximum  p  which  aatiaflea  the  foUcwiag  Inequality 


^  (p-2) 


♦  e' 


1  <0 


1  -  e 


1  >  e 


<3.  32) 


Then  p^^yjh)  are  calculated  aa  a  function  of  T  and'h.  In  Appendix  C»  the 
proof  la  given  to  ahow  that  for  an  unaaturated  eacillatlen  la  alwaya 
equal  to  or  aaaaller  than  _ of  aaturated  eacUlation. 

#IIn( 


Thla  relation  will  be  very  uaeful  for  ttie  analyoia  of  enaaturated  oaclllatlon 
aa  well  aa  for  the  derivation  of  the  gain  boundary  of  auch  oaciUationa.  The 
gala  boundary  of  the  unaaturated  oacillatien  la  given  by  eK^fh)  of  Eq.  (3.  22) 
when  p  takea  the  value  of  p^^^(h)  obtained  above.  Thna 


aK,(p.h) 


J^phlUe 


bph<Ue“*‘®^)+2e 


(e*-l)#|p. 


>bT) 


(3.  34) 


udiere  p  ■  p^^^(h)  obtained  from  Eq.  (3.  f  2). 

The  valuea  of  a  K^(h)  are  calculated  for  variena  h  and  T  and  ahewh  in 
Fig.  (  8).  It  ia  noticed  that  eK^(h)  ia  vary  ctoae  to  cK^fT)  or  alifhtly 
above  aK,(T),  aa  aeon  from  the  aame  figure. 
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FIG. 8  0(Ks(ALmax)  FOR  SATURATED  AND  UNSATURATEO 
OSCILLATION  AND  EXPERIMENTAL  STABILITY 
BY  I.B.M.  704 
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It  can  b«  observed  by  comparing  Eq.  (B.  1)  and  <B.  7)  of  Appendix  B 
that  aK^(p.  T)  le  smaller  than  eK^(p,h)  as  long  as  X)(p)  is  positive. 

aK,(p.T)  <  aK^(ph) 


if  D(p)  >  0 


(3.  35) 


Then  using  the  relation  of  Eq.  (C.  1)  that  has  been  proved  in  Appendix  C.  the 
following  relation  can  be  derived. 

•*^s<**max<^>*  T)<aK>^,Jh),  h)  (3.36) 


for  D(p)  >  0 

D(p)  is  positive  for  bT  <  1. 9. 

Therefore  aK,(p^^^(T).  T)  is  the  lowest  gain  boundary  for  the  existence  of 
saturated  and  unsaturated  oscillation. 

When  bT  >  1.  9.  0(p)  <  0 

hence, 

aK,(p.  T)  >  oK^(K.h) 


And  the  smallest  of  e  K^(|i,h)  is  obtained  by  letting  h  e  0  in  Eq.  (3.  34). 
Eventually  it  is  reduced  to  2b.  Therefore  a  K  ■  tb  yields  the  lowest  gain 
boundary  for  bT  >  1.  9.  Combination  of  these  two  boundary  curves  presents 
the  sufficient  condition  for  nra-existence  of  saturated  and  unsaturated 
oscillation  of  PWM  systems.  This  is  plotted  in  Fig.  (  8)  and  we  can  observe 
that  it  is  very  close  to  the  experimental  data  obtained  by  IBM  704.  We 
mentioned  that  the  pulse  width  of  square  wave  approximation  can  be  regarded 
as  the  average  value  of  the  pulse  width' during  one  period  of  limit  cycle. 
However,  andther  approach  to  this  probUm  is  possible  and  it  is  exptained 
in  the  following.  When  we  find  the  maximum  pulse*width  h  and 

fllWK 

pulse-width  h^^^  during  one  period  of  liaalt  cycle,  we  can  calculate  tts 
corresponding  aK,(h^^)  and 

The  actual  gain  boundary  corresponding  to  tte  actual  wave  shape 
must  lie  in-between  these  two  gains.  But,  we  knew  ^at  is  an  increasing 
fuactien  of  h  and  a  K,||,  h)  is  an  decreasing  function  of  p  as  well  as  of  h. 
Therefore 
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and 


‘‘max  ^**max^  ^  ‘‘max^*‘mln^ 

aK,(»ij.  h)  <  oK^iiij,  h) 

for  (ij  >  |i2  <^<1  >  0 

Combining  these  properties  we  may  conclude  that  aK^  cerrespoadiag  to  h^^ 

is  smaller  than  the  one  corresponding  to  h^j^  for  bT  <  1.  9.  Hence  we  take 

a  K  (h  )  as  the  lowest  gain  for  the  existence  of  the  uneatnrated  eecillatlon 
s'  max'  • 

of  which  maximum  pulse-width  Is  h^^.  Actual  gain  must  he  above  this 
value.  When  bT  >  1.  9.  we  Uke  oK  *  2b  as  the  lowest  gain. 

3.  4  Extension  of  Theories  to  Other  Types  of  Nonlinearities. 

The  noticeable  result  which  has  been  derived  in  this  chapter  is  that 

the  method  to  obtain  u  and  the  corresponding  oK.  can  be  extended  to 

max  s 

other  types  of  nonlinearities. 

Whether  the  nonlinearity  is  a  saturating  amplifier  with  linear  regien 
or  a  quantised  level  amplifier  does  not  matter,  provided  that  they  have  a 
complete  saturation  as  shown  ia  Fig.  3.  We  can  extend  our  theorise  to  the 
analysis  of  such  systems 

The  limitation  on  the  longest  period  of  limit  cycles  for  such  aonlinear 
systems  is  exactly  identical  with  the  one  for  the  PWM  system,  because  they 
will  behave  Just  like  a  relay  system  when  the  esciUation  remains  in  the 
completely  saturated  region.  Therefore,  the  equation  to  find  given  by 

Eq.  (3.  8)  and  Eq.  (3. 14),  is  valid  without  any  change  for  the  nonlinearities 
which  have  the  completely  saturated  region. 

We  need  a  slight  modification  on  the  equation  that  gives  the  gain 
boundary  of  saturated  mode  oscillation. 

Let  a  be  the  tangent  of  the  inclination  cf  the  line  which  combines  the 
origin  and  the  edge  of  the  saturated  region.  Then,  as  far  as  the  saturated 
region  is  concerned,  we  have  the  following  input*outpttt  relations  of  aonlinear 
component, 


(3.  37) 


e;(t)  »  I  l<a«„<0) 

ej^d)  ■  N(e^(0)l  l£-l<a  eJO)  <  I 

■  -1  if  a  e  (0)  <  -1 

B 

N  [  ]  in  the  above  aquation  indicate  a  the  particular  noaliaear  amplitude- 
depending-function  which  can  be  apecified  for  each  caae  of  noaliaearitica. 

In  case  of  auch  amplitude-dependent  nonlinear  functiona,  £q.  (12.  5 ) 
or  (2.^  )  can  atill  be  uaed  to  obtain  the  reaponaea,  and  the  output  of  the 
nonlinear  component  E^^^Ca)  ia  given  aa  followa. 


»  N(eJO)) 


1  -  e 


-Ta 


(3.  38) 


Subatituting  thia  into  Eq.  (2.  68)  yielda 
1 

■ 


^ie 


'T  I  777^  f.„ 


dp 


(3.  39) 


When  the  oacillation  ia  in  the  completely  aaturated  mode,  the  nonlinear  gain 
factor  given  by  N(  )  will  become  equal  to  1  or  •!  aa  ahown  in  Eq.  (3.  37). 
Hence,  the  equation  which  yielda  the  reaponaea  at  the  aampling  iaatanta 
becomea  identical  with  the  one  for  the  PWM  ayatem  given  by  Eq.  (3. 1). 
However,  becauae  of  the  aaturation,  the  input  to  the  nonlinear  component 
during  the  firat  half  period  of  the  limit  cycle  becomea  aa  followa. 


Subatituting  Eq.  (3.  39)  into  the  above  equation  and  obtaining  ita  aum  for  it 
from  o  to  p-1  yielda 


dp 


(3.  «) 


39 


(3.  42) 


Hence  the  fain  boundary  for  eaturated  mode  oacillatioa  becomea 


r 


It  is  easily  observed  by  comparing  the  above  aquation  with  Sq.  (3.19)  lor 
PWM  systems  that  a  for  amplituderdspsadsnt  nonlinoarltiss  can  bo 
obtained  by  multiplying  1/T  to  oK^  for  the  PWM  system. 

This  multiplied  factor  1/T  gives  an  inclining  character istic  to  oK^  as 
shown  in  Fig.  (0  ),  different  from  the  flat  characteristic  of  oK^  lor  the 
PWM  system  which  is  shown  in  Fig.  (  7 ).  a  lor  the  saturating 

gain  amplifier  as  well  as  for  the  quantised  level  amplifier  are  plotted  in 
Fig.  (10)  and  (11).  They  give  sufficiently  close  gain  boundaries  to  the 
experimental  data  on  IBM  704. 
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FIG- II  o^K5(^=I),o^Kjj,AND  mox^  FOR 

QUANTIZED  LEVEL  AMPLIFIER  AND 
STABILITY  BOUNDARY  BY  I  B  M  704 


CHAPTER  IV 


Stability  of  tho  Equilibrium  Point 

It  wao  ihowa  ia  roforoaco  8  that  tho  stability  of  ths  limit  eyclo  is 
tostsd  by  oaamiaiag  ths  sigsavaluss  of  ths  charactsristic  matrix,  which 
coatists  of  ths  first  partial  dsrivativss  of  ths  rariabtss.  Also  ths  squUibrium 
point  can  bs  rsgardsd  as  ths  limit  cycls  of  oas  sampling  psriod. 

Ths  basic  stsps  to  tsst  ths  stability  of  limit  cyclss  which  ars  giroa 
ia  rsfsrsacs  8  wsrs  prsssatsd  ia  Chapter  I.  The  basic  assumption 
ol  this  liasariaatioa  of  ths  aoaliasar  diffsrsacs  squatioa  is  that  ths  psrtur- 
batioa  around  ths  limit  cycls  or  ths  squUibrium  poiat  is  suflicisatly  smaU. 
Howsvsr,  ths  actual  limitatioa  oa  ths  magaituds  of  ths  psrturbatien  uMch 
psrmits  such  appsoxUnatioa  is  aot  spscifisd  sxplicitly.  Ws  wUl  indicate 
that  suck  approximation  bscomss  valid  by  svaluatiag  ths  error  which  is 
brought  about  by  linear  approximation. 

Ws  will  derive  ths  total  liasarisation  tschaiqus  for  UmU  purpose  ia 
general  form.  It  is  known  that  ths  diffsrsacs  equations  of  ths  rsspoass  and 
its  derivatives  for  linear  samplsd*data  systems  ars  written  in  ths  following 
matrix  form. 

'iii  1  •  la’V  1  ''n  1 1  o'i’, )  J  • ».  1 . »•>  (  4.1  ) 

where  consists  of  system  constants  and  rsprsssats  ths  effect  of 
the  input,  and  q  is  ths  order  of  ths  linear  plant.  SimUar  rsprsosatatioa  is 
possible  for  aoaliasar  samplsd«data  oystsms,  using  a  vector,  X_,  for  the 
system  variables  and  a  nonlinear  vector  N(X^^  for  ths  output  from  ths  nonlinear 
component. 

•  10)  ♦  W*,!!  (4.2) 

and  it  is  assumed  that  X^  and  Vl  is  on  ths  squUibrium  poiat,  8^,  which 
is  ths  origin  ia  most  cases.  Then  if  we  lot  tho  compsaoats  of  ths  vector, 

X^,  bs  Xjj^,  x^^,  . . . ,  x^^,  these  x^^^  corresponds  to  c^^  ia  Bq.  (  4*  I  )• 

Also  ths  charactsristic  ssatrix  iQ]  of  ths  liaoar  plaat  corresponds  to 

U.l  -)- 

Letting  with  its  compeasat  y^^,  . . . ,  y^  bo  ths  small 
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disturbance  of  X  ,  Eq.  ( 4.  2  )  is  rewritten  when  such  disturbance  is 
n  n 

added  to  X  . 

n 

♦  VI  *  *  ’'n'l  '  ■>.  3  > 

or 

Vl)"  Vll*  101  V)  V)  > 

However  in  most  nonlinearities,  the  output  of  the  nonlinear  component  is 
only  the  function  of  input  itself,  and  is  not  the  function  of  its  derivatives. 
Therefore, 

N(X^)  =  N{XjJ  >  N(-c^)  (4.5) 

In  that  case,  all  the  partial  derivatives  of  N(X^)  vanish  when  it  is  differentiated 
by  for  j  ^  1.  Hence  N(X^  Y^)  can  be  expanded  into  powers  of  x^^ 

by  Taylor  expansion 

*  ^n^  “  V  ^  ^In  ^  TT  7^  +  •  •  •  (  4.  6  ) 

-in  axj^ 

where  all  the  partial  derivatives  are  evaluated  at  X^  ■  S^.  Denoting  the 
summation  of  all  the  terms  as  N'(y|^)  except  for  the  first  term  in  the  above 
equation,  it  becomes 

Substituting  this  into  Eq.  (  4.4)  and  using  Eq.  (  4.  Z  )  yields  the  difference 
equation  of  the  perturbation  around  the  equilibrium  point. 

Vil  ■  loiV^N-ly,,!)  u.e) 

Then  denoting  the  summation  of  all  the  terms  as  N"(y^^)  except  for  the 
first  two  terms  in  Eq.  (4.  6.1),  we  have 

V  ♦  ""'W  '  ’> 
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Substituting  this  into  Eq.  (  4.  8)  yislds 


IQI  \]  + 


-  lAj  t  J  4  N»(yj,)J 


(4.10) 


whsrt  ths  compoasat  ia  the  first  column  of  is  ths  sum  of  and 

the  corresponding  term  ia  Q.  and  tho  other  components  of  tiro  matrices  are 
all  identical. 


a 


il 


^ij  *  J  •  ,  q  ( 4. 11 ) 

a 

It  is  observed  that  ( A^]  given  above  is  identical  with  I A^]  of  Eq.  (144)  in 
the  previous  section.  However,  Eq.  (  4. 10)  has  no  appreadmation  eriiereas 
Eq.  (i.l4  is  the  first  approximation  of  Eq.  (^.4fi),  aeglectiag  aonlinear 
vector,  N"(y|^)].  For  certain  types  of  aoalinearities  it  can  be  proved 
that  N"(Y.)  is  finite  inside  certain  region  S,  which  is  a  vector  space  of  Y.. 
In  other  words,  we  can  find  such  and  Rfwith  thsir  components 
and  R^,  respectively)  that 

N'i<yin»  ■  ^^n^ln  «or  aU  i  (4.12) 


end 

0  <  R|^  <  R^  for  all  a  and  i  if  Y^e  8  (  4.'t3) 

R^<io  a  positive  finite  number  and  so  is  R^  .  The  number  h  with  its  component 
is  introduced  for  mathematical  convenience  to  represent  the  sign  and 
gain  factor  and  also  is  a  finite  number.  Then  substituting  Eq.  (  A  12)  into 
Bq-  (4»io)yi«w* 


Y 


lA.)  \]  ♦  k  R, 


<4.  lit 


where  the  compoaeal  in  the  first  column  ol  is  the  aum  of  the  corrospanding 
component  of  A^  and  the  component  of  kR^,  and  the  oihor  companonts  of 
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two  matrices  are  all  identical. 


*il  *’  ^i*^in 


Qil  ^4^/  ^  *^in 


‘"ij  *  ‘ij  “  % 


for  j  «  2 . . 


U.15) 


This  linearised  difference  equation  in  its  form  shown  in  Eq.  (4.14*)  i* 
essentially  a  nonlinear  difference  equation  because  the  matrix  [B^]  is  a 
nonlinear  matrix  containing  the  nonlinear  factor,  R^.  However,  Eq.  (  4,14 ) 
is  derived  without  introducing  any  approximation  and  is  valid  inside  the 
vector  space  S  for  all  n.  And  we  know  that  R  is  bounded  as  seen  in  Eq. 

(  4  13)  for  all  n,  hence,  applying  the  stability  criterion  to  the  stxatrix 
[B  ]  in  the  same  manner  for  the  lineag  matrix,  we  may  derive  certain 
conditions  fox  the  stability  of  the  difference  equation  of -Eq,  (4.14!),  using 
the  boundedness  of  R  .  When  the  condition  obtained  above  i#  satisfied, 
all  the  eigenvalues  of  lie  inside  the  unit  circle  for  all  n. 

Although  this  condition  does  not  necessarily  guarantee  that  any 
disturbance  — a  0  as  n  —a  00,  it  gives  stricter  restriction  on  the  stable 

region  than  the  condition  on  the  incremenally  linearised  aaatrix 
which  is  derived  in  the  previous  section.  Also,  by  evaluating  explicitely 
the  magnitude  of  the  error  factor  kR,  which  is  caused  by  the  liaesxisatien, 
we  can  estimate  the  range  of  S  in  which  the  totally  linsariasd  matrix 
[B^]  can  be  regarded  as  a  time* invariant  matrix  [B].  Within  that  region, 
the  test  on  the  single  matrix  [B]  ensures  the  stability  of  the  equilibrium 
point.  Then,  if  [B^]  is  a  stable  matrix  for  all  a,  i.  e. ,  if  the  eigenvalues 
of  (B^J  lie  inside  the  unit  circle,  the  equilibrium  point  to  stabts  within 
the  region  that  the  matrix  [B^]  can  be  regarded  as  a  constant  matrix  (B], 
Such  a  region  can  be  found  by  evaluating  the  magnitude  of  the  error  vector 
kR  and  by  comparing  it  with  other  time  invariant  components  of  the  matrix 

®n- 

Example  4. 1 

The  total  lineariaation  method  is  illustrated  on  the  PWM  system  and 
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its  ussfulnsss  will  bs  provsd.  The  linsar  plant  chosaa  art  again  first- 
and  second'ordar  systems. 


(1)  First-Order  System.  For  the  first-order  system,  the 
difference  equation  is 


where 


h  «  -oy(a)  c. 


(  4. 16) 
(4.17) 


Obviously  X_  ■  x,  ■  c  and  Q  ■  and 

n  lA  n 


-bT  ,.**n 


N(XJ  .  N(c^)  .  |,{n)e-**  (e  “-!) 


(4.1# 


in  Eq.  (4.  The  equilibrium  point  ■  0.  The  difference  equation  of 
the  perturbation  y^  around  the  origin  becomes,  from  £q,  (  4.  JQ)  and  (  4.  n) 


'ntl  *  Vn*  ""‘I'.l 


where 


and 


A 


i-2  n 


c  *0 
n 


_*bT  if.'bT 
»  e  -uKe 


(4,17^) 


(4,j5« 


00 

-  ^  7T-(-abir(n)  y^)*  (4,  a) 

iwZ 


By  Eq.  (  4.12) 


N”(y„)  .  kR^y„ 

where  k  is  chosen  as  follows 
k  >  -aX 

Then  equating  Eqs.  (4.  21 )  and  <4.  22).  we  find  R 


(4.  22) 
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(  4.  22> 


1*2 

U  U  notie«d  that  -^(a)  !•  alwaya  positlva  bacaaaa  y<a)  ■  -1  tdwa  y^  >  0 

•«!  Y<«)-  •fl  whaa  y^  <  0.  Alao  iay^|  will  aot  axeaad  tha  aam^lag  pariod 
T  bacauaa  af  aaturatloa.  Thua 

0  <  -  a  y<a)  y^^b  <  bT  (‘4.2<^ 

Thaa,  uaiag  Eq.  (4.  23)  aad  tha  aboaa  reault,  ia  booadad  aa  Ifllowa: 

ia2 

■  CT  1  ^5) 

-  R 

Tha  tatally  liaaariaad  dilfaraaca  aquation  of  tha  parturbatioa  arouad  tha 
origia  ia  from  Kq.  (4. 14) 

W  •  V.  M  »» 

whara  la  obtaiaad  from  Eqa.  (4.15)  aad  (4.  20) 

n 

m  ♦  kR^  +  (1  .  oKla"^"^  .  oKR^  (4,  27) 

Eq.  (4.  26)  ia  ataUa  if 

1  B,l  «  I  (4.  Ml 

•uhatitutlag  Eq.  (4. 27)  lato  tha  abova 

|a’^’^(l.aK)-oKR^|  <1  (4.29) 

Thia  ia  rowrittaa  aa  lollowa: 

(I  -  a’^^  ♦  a  Ka’^'^  ♦  o  KR^Kl  ♦  •  a  ♦  R^)  )  >  0  (4.  2f) 
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The  content  of  the  first  brocket  is  positive  since  oK  and  R  are  aon«negative 

B 

quantities.  Hence,  the  stability  condition  is  reduced  to 
•bT 

a  K  <  ^  (A.  31 ') 

is  not  a  constant,  but  is  restricted  between  0  and  R  for  all  a  by  Eq.  ( 4.25)> 
Therefore,  the  maximum  value  of  R.  is  substituted  into  Eq.  (<4.31>)  ia  order 
to  obtain  the  lowest  boundary  of  oK.  Thus 

•bT  •bT 

aK  <  ^  ^  ”.'bT- 

R  >  e  I  -  e 

Various  important  results  have  been  obtained  concerning  the  stability  of 
the  first-order  PWM  system.  Reviewing  and  combining  thedb  results,  a 
conclusion  of  the  stability  of  ouch  system  can  be  derived. 

In  the  first  place,  it  is  proved  that  the  longest  period  of  relay  mode 
oscillation  is  two  sampling  periods,  p«l.  In  the  second  place,  the  gain 
boundary  corresponding  to  p  s  1  is  obtained. 

oK,(p  s  1)  «  bT 

Thus,  we  may  state  that  no  relay  mode  oscillation  can  exist  below  aK^(p  ■  1). 
In  the  third  place,  it  has  been  proved  that  the  origin  is  stable  for  any  pertur¬ 
bation  below  the  gain  which  is  given  above.  This  fact  eliminates  the  possi¬ 
bility  of  the  existence  of  PWM  oscillations  or  any  other  irregular  oscillations 
below  o  K^(p  s  1). 

By  these  three  steps,  we  could  have  successfully  eliminated  the 
existence  of  all  the  types  of  oscillations  below  ■  1).  On  the  other 

hand,  it  has  been  shown  in  the  previous  section  that  the  relay  mode  oscillatien 
can  exist  above  aK^(p  s  1)  and  ie  stable.  Hence,  we  have  shown  that 
oK  <  oK^(p  «  1)  is  the  necessary  and  sufficient  condition  for  the  absolute 
stability  of  the  first-order  PWM  system. 
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(2)  Sccond«Ord«r  SytUm 

Next  the  case  of  the  second* order  system  will  be  discussed.  The 
transfer  (unction  is  given  by  Eq.  as 

and  the  difference  equations  will  be  derived  M  iadliwra 


„  -b<T-h  )  . «  , 

'n.l  •  V'"l  ^  '  “n  •  •  ‘  '-4 

- 

n  n 


-bT 


(4.  U) 


K  . -bT 


e  “*) ♦  e 


•bT 


«<"n‘ 


(4.  34) 

Then  the  components  matrix  ( A^l  of  Bq.  pUdl  is  givsa  by  Eq.  U.15). 

l^ol  • 


*11  *12 


‘21 


‘22 


(4.  35) 


where 


II 


12 


a„  a 


*21 


22 


•f 

<  4.  3d 

sr 

1 . 

(4.37') 

• 

“-B - 

-NT-hn) 

K' 

-  oKe 

(4.38) 

•JL 

.-bT 

•  • 

(4.39) 
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Let  =  c^,  X2^  =  «nd  be  the  corresponding  perturbations, 

and  the  equilibriunn  point,  S^,  be  the  origin.  Then  [Q]  X^]  of  £q.  (4.  2.  ) 
represents  first  two  terms  in  the  right  hand  side  of  Eqa.  (4.  33)  and  (4.  34) 
and  N(X^)  represents  the  last  terms  of  two  equations.  Obviously  N  is 
solely  a  function  of  c^,  since  h^  s  -a  ^(n)c^.  The  components  of  Eq. 
(#,  10  )  are  given  by  Eqs.  (4.  35)  to  (4.  38).  Choosing  ■  oK,  k2  ■  -oK 
in  Eq.  (4. 12)  we  can  find  and  R^n' 


N'ilW  •  ^  IT 


oo 

S  4. 

,(i) 

L  TT 

is  2 

n 

'I  i 

—  yin 


(4.  40) 


c  s  0 
n 


and 


oo 

i=2 


8c 


nr"  Mn 


(4.41) 


c  =  0 
n 


Then  the  limiting  value  of  Rj^  and  R2j^  are  obtained. 
From  Eqs.  (  4.  3)  and  (4.40)  we  have 


aKR 


aKe 


-bT 


In 


oo 

^  (-aby(n)yj^) 


i-l 


Ix. 


Because  of  the  saturation. 

0  <  -Q  b'Y{n)yj^  <  bT 

Hence 


(4.  42) 

(4.43) 


(4.44.) 
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Rj  =  R/b 


where  R  ie  given  by  Eq.  (4.  21).  Similarly  from  Eqe.  (4.34)  and  (4.41  ) 

00 


.  -a 


i«2 


(4.45) 


Hence,  referring  to  Eq.  (4.44) 


°  -*2n  ■  **^n 

The  totally  linearised  form  is  given  by  Eq.  (  4.14') 

Vl  1  •  t®n> ’'.1 

where  the  components  of  are  given  by  Eq.  (4. IS)  and  Eqs.  (4.  35)to  (4.3t) 


‘12' 

m 

•>12  ■ 

•21  ■  •“*20  . 

*22. 

“21 

•>22 

C_  »  0 

b 

(4.48) 

is  stable  if  the  following  two  conditions  are  satisfied  by  its  components?^ 
L  (I  +  bjibjj  -  bijbji)(l  -  biibjj  +  bjjbji)  >  0  (4.49) 

n.  (1  +  bjjbjj  -  bjjbjj  +  bjj  +  *>22X1  +  ^11*^22  '  **12**2l  *  ^11  ’  **22^  ^  ® 

(4.50) 

Condition  I  is  tested  first,  by  using  Eqs.  (4.35 1  to  (4.38  )  and  (4.48) 
^11^22  ■  **12^21  ■  •  « *^n  (4.  51  ) 

The  content  of  the  first  bracket  is  clearly  positive  sines  is  a  non^nsgative 
quantity  by  Eq.  (4.44).  Hence,  the  content  of  the  second  bracket  must  be 
positive  to  satisfy  Condition  L 

1  >  e***’’  +  oKRj^ 
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aK  < 


1  -  • 


(4.52) 


The  above  inequality  ie  eatiafied  for  all  a.  if  i»  replaced  by  its  maximum 
value  which  is  given  by  Eq.  (4.44) 


oK  < 


1  -  e 


(1  - 

1  -  bTe'®'^  -  e 


(4.53) 


Thus  Condition  1  is  satisfied  for  all  a  K  <  a  K . 


(4.48) 


Next,  Condition  II  is  investigated.  Using  Eqs.  (4.35)  to  ( 4.38)  Md 


^  *  **11^22  ‘  **12^21  “  *  **22^ 


(4.54  ) 


and,  also 

^  ^  ^l‘*22  ’  ^2^21  *  hi  *  **22 

I 

«  2(1  +  e'***^)  -  (1  -  e’**’’)  +  ZoKRj^ 

and  this  must  be  positive  to  satisfy  Condition  II.  Hence, 


(4.  55i) 


a  K  (  1-1^ -  -  2R2^)  <  2(1  +  e"***^) 


(4.56) 


This  inequality  is  satisfied  for  all  n  when  R2^  is  replaced  by  its  minimum 
value  which  is  eqqal  to  sero  by  Eq.  (4. 46).  Thus, 


oK  <  2b 


(1  e 
1  -  e’ 


■  a  K. 


n(4.57  ) 


It  is  observed  that  this  o  is  identically  equal  to  the  lowest  gain  for  the 
PWM  oscillation  of  M  ■  2  which  is  obtained  in.rjtference  X  When  these  two 
conditions  are  combined  we  may  state  that  the  origin  is  stable  if 


( 


qK  <  min(aK^.  aK.) 

C  Q 


and  if  the  system  remains  within  the  region  that  the  matrix  [B^]  can  be 
approximated  by  a  tin.e  invariant  matrix  [B]. 

Thus  we  could  have  succeeded  in  eliminating  all  the  regular  or 
irregular  oscillations  around  the  origin. 

In  the  first  order  system,  the  test  on  [B^]  gave  the  absolute  stability 
boundary  oK^. 

In  the  second  order  system,  we  could  have  derived  a  from  Condition 
I  which  could  not  have  been  obtained  from  the  incrementally  lineariaed 
matrix  [A  1. 

^e  have  mainly  discussed  the  saturated  oscillations  in  Chapter  XII, 
and  the  small  oscillation  has  been  studied  in  this  chapter. 

Combining  the  results  obtained  from  these  two  different  approaches  will 
give  a  good  means  of  solution  for  the  stability  problem. 

and  the  experimental  stability  curve  obtained  by  IBM  704 
are  shown  in  Fig.  (8  ).  Also  min(aK^,  oK^)  are  plotted  in  Fig.  (1,9) 
together  with  the  experimental  curve.  It  can  be  observed  that  both  curves 
give  sufficiently  close  values  to  the  experimental  data. 

The  stability  condition  obtained  by  Lyapunov's  method^  is  added  to 
these  figures  for  the  sake  of  comparison. 
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FIG.  9  Q(K 


aXd  and  experimental  stabil itv  o«.. 

BV  IBM  704  ^  boundary 


CHAPTER  V 
CONCLUSION 


5.1  Conclusion 

Periodical  oscillations  within  nonlinear  sampled>data  systems, 
as  a  whole,  has  been  investigated  in  this  report. 

The  fundamental  equation  that  gives  the  exact  feature  of  limit  cycles 
will  be  very  useful  for  the  precise  analysis  of  such  systems.  Oscillations 
which  are  almost  periodical  can  also  be  treated  by  this  equation  for  an  ap¬ 
proximate  description  of  their  behavior. 

Thus  the  fundamental  equation  and  its  modified  equation  can  be  very 
powerful  tools  in  investigating  the  steady  state  of  nonlinear  eampled-data 
systems.  When  the  periodical  input  is  applied,  the  same  equation  yields  the 
desired  responses,  and  the  results  are  superior  in  their  precision  to  the  one 
obtained  by  the  describing -function  method.  They  are  based  on  the  transform 
method,  and  can  be  applied  for  any  order  of  the  plant  transfer-function.  More¬ 
over,  they  can  be  used  for  any  sl4ape  of  nonlinearitiee,  not  necessarily  PWM, 
but  also  relay,  saturating  amplifier  and  quantised  level  amplifier. 

Similar  approaches  to  the  problem  (f  limit  cycles  are  possible  as 
shown  by  Shao  Oa  Chuan  ^  by  means  of  the  canonical-form  represeatation  or 
by  H.C.  Tomg  means  of  the  discrete-function  method.  However,  in  these 
studies,  setting  up  of  state  equations  or  difference  equations  for  each  case  is 
always  required  and  the  siae  of  system  equations  will  become  larger  and 
larger  when  the  order  of  the  plant  increases. 

On  the  other  hand,  in  our  transform  method,  every  equation  that  gives  the 
exact  feature  of  responses  can  be  derived  from  one  fundamental  equation 
which  is  common  to  all  the  plants  and  to  all  the  nonlinearitiee. 

The  steps  that  must  be  followed  are  Just  substitutiag  the  actual  form  of 
the  plant  transfer -function  into  the  fundamental  equation  and  giving  considera¬ 
tion  to  the  particular  type  of  nonlinearitiee. 

The  stability  problem  was  another  important  topic  of  this  thesis.  As 
was  mentioned  in  Chapter  I,  our  attitude  towards  this  problem  is  microscopic 
in  contrast  to  the  macroscopic  approach  such  as  Lyapunov's  second  method. 
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Pimiafctiaf  th«  pottlbU  •xitUnca  of  oil  Mturotod  mod  oaiotarotod 
oicUlotlOBO,  tho  ■olAciont  condition  for  tho  noa-oxiotoaco  of  ouch  otcillatloat 
la  tho  PWM  iyotom  ha*  boon  dorlvod  and  lo  eomparod  with  tho  oaporimoatal 
roottlto.  ao  wall  as  with  tho  oufficioat  condition  for  tho  asymptotic  stability 
in  tho  largo  as  dorlvod  by  Lyapunov's  socond  mothod. 

Tho  sti^ility  boundary  derived  by  our  mothod  givos  a  closer  criterion 
for  tho  oxporlmeatal  results  than  that  obtained  by  tho  Lyapunov  mothod.  The 
roaoon  for  this  can  be  attributed,  in  the  first  place,  to  the  difficulty  in  finding 
the  boot  Lyapunov  function  for  the  specified  types  of  nonlinearity;  in  the  second 
ptaco.  lo  the  fundamental  property  of  auch  a  macroscopic  method  wherein  ono 
hwot  oapoct  that  the  worst  case  might  happen,  without  paying  attention  to  tho 
liasliatlon  on  oscillations  inside  the  system  imposed  by  tho  operation  of  the 
nonlinear  function. 

We  would  never  deprecate  tho  approach  from  the  macroscopic  point 
of  view.  Wo  might  oven  expect  that  tho  straightforward  method  to  find  tho 
boot  Lyapunov  function  to  give  the  necessary  and  sufficient  condition  for  tho 
stability  could  bo  estaUlshod.  However,  we  believe  that  out  localised 
approach  is  also  useful  in  do  signing  nonlinear  sampled*  data  systems. 

It  is  frequent  that  tho  system  designer  wishes  to  eliminate  tho  partleular 
aaodoo  of  oscillations,  especially  of  tho  fundamental  frequency  (half  of  tho 
aaaapliag  frequency)  or  of  a  few  of  its  subharmonics.  Or  ho  may  attempt  to 
iliasiaali  only  tho  oocillatiens  of  largo  amplitude  which  probably  remain  in  tho 
oalatatod  regions,  hi  such  cases*  tho  mothod  that  wo  have  ostabliohod  is 
directly  applicable;  and  the  designer  will  obtain  the  sufficient  and  satisfactory 
iafomsatioa  on  his  problem  without  wasting  time  in  the  struggle  to  find  a 
eamprohoaoivo  stability  condition  that  may  lead  to  an  rxcessively  conservative 
rooult. 

Tho  modification  of  Noaco'smothod  to  tost  tho  stability  of  an 
Ofuilibrium  point  is  attempted.  It  is  worthwhile  to  notice  that  tho  stability 
boundary  obtained  by  this  mothod  yields  tho  closed  cr  itorlon  to  tho  oxpori- 
nMuml  result  as  shown  in  Fig.  9 

Bnsphaols  is  placed  on  tho  PWM  system  in  this  thesis,  although  the 
aimlyeiB  has  alwayo  boon  oxtondod  to  other  typos  of  nenlinoaritioo. 

As. Air  as  tho  analysis  is  concornod,  ao  approximation  is  introduced* 
ouah  ao  small*  oipMd  condition  or  oino-wavo  approximation.  Tho  PWM 
oontroUor  is  a  vary  sonoitlvo  device  while  it  is  operating  in  tho  unoaturatod 
region)  also  by  moans  of  Its  saturating  property  it  prsvonts  tho  plant  from 
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r«c«lviag  an  •xc«*iiva  influanca  from  tha  input. 

Wa  obaarve  in  Fig.  9  that  tha  oparatlng  ragloa  (atabla  ragioa)  of 
tha  PWM  ayatama  la  rathar  uailorm,  daa  to  tha  flat  eharaetarlatlc  of  tha 
atabllity  boundary,  Thar af ora.  tha  ranga  of  tha  choica  of  aampliag  frafuanciaa 
ia  vary  wlda.  On  tha  othar  band,  in  caaa  of  othar  aoalinaaritiaa,  tha  choica 
of  aampliag  fraquancy  ia  rathar  limitad  bacauae  of  tha  Incliaad  charactariatica 
of  tha  atabUity  boundary  (Figa.  10,11.  Tha  PWhf  ayatam  ia  auparlor  to  tha 
ralay  ayatam  in  tha  aanae  that  tha  formar  haa  a  atabla  ragioa  whila  tha  lattar 
can  never  get  rid  of  the  oaciilationa  of  the  fundamental  frequency  and  ita 
firat  aubharmonica.  Alao  it  la  obaarvad  axparimantally  that  tha  PWM  ayatam 
reachaa  the  equilibrium  atate  conaidarably  faatar  than  tha  othar  typoa  of 
aoalinaaritiaa  auch  aa  the  aaturating  amplifier  with  linear  region,  otarting 
from  tha  aama  initial  conditimia.  From  thia  fact,  wa  nay  atata  that  tha 
PWld  cantroUar  ia  a  more  aan.aitiva  device  than  othar  nonlinear.  aiApiitttda- 
dapandant  controlling  davicaa.  And  thia  will  offer  an  advantage  to  tha  optintal 
control  by  maana  of  tha  PWU  ayatam. 

S.  2  Sttggaationa  for  Future  Work 

When  tha  nonlinear  oamplad*data  ayatam  ia  rapraaantad  by  a  aat  of 
nonlinaar  diffaranca  aquationa.  it  noay  be  reduced  to  a  totally  linaariaad  form 
aa  followa: 

Y(n+1)].  [B„)Y(n)] 

where  Yfa^  ia  a  vector  rapraaanting  tha  atata  at  tha  nth  aampling  inatant  and 
[Bq]  ia  a  totally  linaariaad  matrix  and  ia  a  nonlinaar  function  of  Y(n^ 

Wa  may  find  a  certain  region  of  B„  where  Bg  ia  atabla.  i.a. .  tha 
aiganvalua  of  [Bq]  liaa  inaida  the  unit  circle  if  Y(n)]balonga  to  tha  certain 
vector  apace  S. 

23 

It  waa  pointed  out  by  Kodama  that  it  ia  falaa  to  eay  that  tha  equi¬ 
librium  point  ia  aaymptotically  atabla  in  tha  large  if  [Bn  ]  ia  atabla  for  all 
tha  pointa  of  atata  apace.  Tha  above  fact  ia  damonatratad  by  hia  counter 
axaaplaa.  by  ahaping  the  nonlinear  gain  curve  for  that  purpoaa. 

However,  wa  believe  that  under  certain  conditiooa  tha  above  will  hold 
true.  And  it  ia  worth  whila  to  look  forward  to  ouch  conditiona.  Thia 
poaaibility  la  aupportad  by  tha  following  fhet.  Whan  tha  naatricaa  [Bj  ] 
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and  [B2  ]  are  stabla,  and  rB2  ara  alao  atabla  by  tha  Frobaniua 

thaoram.  Howavar,  nothing  ia  aaaurad  on  tha  atability  of  tha  product  matrix 
[Bj  ]*[B2  ]•  But  thaaa  two  matrlcaa  ara  cloaaly  ralatad  to  aach  othar  by  tha 
nonlinaar  ayatam  aquation.  If  thla  ralationahip  ia  praciaaly  invaatifatad.  wa 
may  find  tha  conditiona  undar  which  f  Bj  ]*fB2  ]  becomaa  alao  atabla.  And  it 
may  ba  axtandad  to  tha  caaa  whan  tha  numbar  of  multipliad  matt  icaa  ia 
incraaaadt  poaaibl  to  infinity. 

Alao  wa  noticad  that  varioua  concluaiona  that  ara  darivad  for  individual 
matrlcaa  cannot  ba  appliad  to  thair  product  at  all.  In  ganaral.  wa  hava  to 
carry  ovar  tha  tadloua  multiplication  of  matrix  alamenta  and  wa  uauaUy 
find  a  vary  dlffarant  concluaion  from  what  wa  axpact  from  tha  individual 
matrix.  It  will  be  vary  halpful  if.  for  axampla.  tha  atability  eritaria  on 
tha  individual  matrix  can  ba  axtandad  to  tha  product  of  matrlcaa  without  carrying 
over  tha  actual  nnatrix  multiplication  proceaa. 

Anothar  auggaation  on  tha  atudy  of  nonlinaar  aamplad-data  ayatama  ia 
tha  appropriate  uaa  of  high-apaad  digital  computora.  Nowadaya»  tha  digital 
computer  ia  extan  aively  adopted  aa  a  controlling  device  of  aamplad-data 
ayatama.  However,  it  can  be  uaad  for  aimulatiag  the  aampled-data  ayatam 
in  ita  programa.  Wa  can  perform  any  type  of  experiment  on  aamplad-data 
ayatama  that  are  incorporated  in  programa  of  digital  computora  in  tha  form 
of  diffarance  aquationa. 

Wa  ballava  that  a  complicated  concluaion  obtained  by  certain  theoretical 
inveatigationa  can  be  accepted  if  it  ia  formulated  for  auitable  uaage  of  digital 
computora,  aince  tha  digital  computer  can  give  the  deaired  data  inatantly 
when  it  ia  wanted.  It  will  ba  adviaable  for  reaearchera  in  aamplad-data  ayatama 
alwaya  to  keep  the  poaaibla  utilisation  of  digital  computora  in  mind. 


APPENDIX  A 


EXPERIMENTAL  WORKS  BY  DIGITAL  COMPUTORS 

1m  ordar  to  verify  the  theoretical  works  developed  in  thie  repor)  high¬ 
speed  digital  computers  are  used  extensively. 

Because  of  the  b%aic  property  of  the  sampled-data  systems  that  has  been 
the  main  subject  of  this  thesis,  the  digital  computer  is  well  fitted  for 
eyntheaUiag  the  system  and  performing  the  experimental  works.’ 

The  function  of  the  closed  loop  nonlinear  sampled-data  system  as  shewn 
in  Fig.  U)  i*  completely  represented  by  a  set  of  nonlinear  dilfarenca 
equations  as  Eq.  (4.16)  or  Eqs.  (4.33)  and  (4.34).  These  nosdinear  difference 
eqioations  can  be  easily  incorporated  into  programs  of  digital  computers 
because  of  their  iterative  properties.  Responses  at  the  end  of  every  sampling 
period  can  be  used  as  the  initial  conditions  for  the  following  sampling 
instants;  this  operation  is  conveniently  performed  by  using  the  transfer 
command  and  the  index  register  in  digital  computers. 

Setting  of  initial  conditions  is  quite  arbitrary  and  the  accuracy  of  the 
computation  is  incompatible  with  that  of  analogue  computers. 

A  series  of  experiments  is  performed  using  Bendix  15  to  verify  the 
existence  and  behavior  of  limit  cycles  as  derived  in  Chapter  H.  Typical 

I 

oscillations  of  PWM  mode  or  relay  mode  are  observed  when  suitable  gain 
is  given  to  the  system. 

oK^  or  oK^  are  accepted  as  boundary  gains  and  a  conspicuous  dif¬ 
ference  of  behavior  of  responses  is  observed  on  both  sides  of  these  critical 
gains. 

Finally,  a  series  of  stability  tests  is  performed  in  order  to  endorse 
the  various  stability  coundaries  obtained  in  Chapter  UI  and  Chapter  IV. 

A  wide  range  of  initial  conditions  (40-b0  points  distributed  on  the  phase 
plane)  is  selected,  and  for  every  initial  condition  the  responses  are  cal¬ 
culated  up  to  100-200  sampling  instants  in  order  to  examine  the  coovergonce 
of'respmses. 

Varying  the  gain  ^  gradually  over  the  critical  regions,  the  border  line 
between  stable  and  unstable  regions  is  traced  with  good  accuracy.  Those 
tests  are  repeated  varying  the  sampling  period  T  and  the  results  ars 
plotted  as  a  function  of  the  sampling  frequency  as  shown  in  Ffg.  } 
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Mnd(9)> 

As  ftB  sxampU,  the  flow  chart  of  the  program  for  the  secoo<t  order 
PWM  system  is  shown  in  Fig.  (12). 

This  is  programmed  for  the  purpose  of  obtaining  the  responses  of 
the  second  order  system  at  every  sampling  instant  when  the  initial  eon* 
dltions  are  specified. 
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FIG.  12  FLOW  CHART  OF  PROGRAMMING  ON  BENDIX  15  (2nd  ORDER  PWM  SYSTEM) 


APPENDIX  B 


PROOF  OF  DECREASING  CHARACTERISTIC  OF 


The  objective  of  this  appendix  is  to  prove  the  gain  boundary  oK^  for 
laturated  oacillation  ie  a  monotonically  decreaaing  function  of  p. 
aK  (p)  it  given  by  Eq.  (3.  23)  for  the  eecond-order  ayatem,  it  can  be  rewritten 
aa  foUowa. 


aK,(p) 


1  4-  2(1 


2b 


vn - 


(B.l) 


In  order  to  prove  the  decreaaing  charncteriatic  of  aK^(p),  it  ia  enough  to 
prove  the  increaaing  charncteriatic  of  the  function  inaide  the  large  bracket  of 
the  denominator  of  the  above  equation. 

Let 


D(p) 


■Jibul.bl]' 

hpTd  +  e-^^^ 


) 


(B.  2) 


then,  the  incremental  difference  AD(p)  ia  given  by 

AD(p)  «  D(p+1)  -  D(p)  (B.  3) 

If  AO(p)  ia  poaitive  for  all  p,  0(p)  haa  a  monotonically  increaaing  chnracteristic. 
Subatituting  Eq.  (B.  2)  into  (B.  3)  yielda  AD(p) 


AD(p) 


bpT(l  +  e">*°^)  . 


-bT)  - 

M;;i)Tji;e^^^*»‘»T, 


)f<ii.  -bT) 


(B.  4) 


f(p,  'bT)  ia  obtained  by  aubatituting  -bT  into  the  place  of  pT  in  Eq.  (3.17) 
or  in  Eq.  (3. 18).  \%hen  thia  ia  aubatituted  into  Eq.  (B.  4)  the  numerator  of 
AD(p)  can  be  reduced  to  the  following  form. 

Numerator  of  AD(p)  .  2  ^  (f♦)(e*^‘*’'-e‘^^‘**^***^) 

tmO 

f«0 
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This  is  obviously  positivs  for  all  p  which  is  s  positive  iatogsr. 
Therefore 


>0  for  p  «  1.  2.  3,  ...  (B.6) 

Thus  we  could  prove  the  moaotmic  iacressiag  property  of  Dip),  hence  the 
monotoaic  decreasing  property  of  aK^(p).  A  similar  property  can  be 
proved  for  aK^(p)  of  uneaturated  oscillations.  It  is  given  by  Sq.  (3.  34)  and 
is  rewritten  ae  follows 


eK^(p,h)  « 


2b 


1 .  t  M.  1 

e“  -1  phTIl  +  e^®*) 


(B.  7) 


pbT(l  ♦  e^—) 

The  function  inpide  of  the  large  bracket  of  the  denominator  of  the  above 

equation  is  related  to  D(p)  of  £q.  (B.  2)  as  follows 

bh 

Function  inside  bracket  ■  ^  (B.  I) 

e**  -  1 


Therefore,  when  0(p)  hae  an  increasing  property,  naturally  the  above  function 
possesses  the  same  property,  which  makes  aK^(p,  h)  have  a  dscroasing 
property. 
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APPENDIX  C 


PROOF  OF  THE  RELATION  u  .  <h)  <  u  (T) 


In  order  to  prove  the  relation  that 

^max  for  0  <  h  <  T  (C.l) 

where  ie  the  maximum  positive  integer  that  satisfies  the  inequality 

of  Eq.  (3. 14)  and  is  the  maximum  positive  integer  that  satisfies  the 

inequality  of  Eq.  (3.  33).  We  denote  the  left  hand  side  of  Eq.  (3.  32)  as  g<h) 

g(h)  -  gj(h)  -  g^th)  (C.  2) 

where 

il(f*)  ■  (C.  3) 

Kk  1  t 

gjih)  «  (e  -1)  [  - — ^ -  1  J  (C.  4) 

I«e  1+e  ^ 


If  the  following  relation  is  proved. 


g(h)  >  g(T)  for  h  <  T 


(C.  5) 


then  it  implies  the  relation  of  Eq.  (C.  1)  since  g(h)  is  a  monetonically 
increasing  function  of  p. 

From  the  following  relation 


82(M 

iTTr- 


h 

T  • 


we  obtain 


(C.t) 


-  ijrrr  -  i^rr  ic-  n 

Let 

gj(H) 

gjHT  •  *  ■  “  “  -  °  •) 


p  >  0 


(C.  9) 


f 


% 


g2(h) 

Tim 


*  I  -  p 


Then 

0  <  a  <  p 


Therefore 

gl(h)  -  gj(T)  -  (gjCh)  -  )  •  -  a(gj(T)  -  |  g2(T) )  >  -o (gj(T)  -  g2(T)  )  >0 

(C.IO) 

since 

gj(T)  -  g2(T)  .  g(T)  <0  for  |i  < 

Eq.  (C.  10)  implies  £q.  (C.  5).  thus  the  proof  has  bssa  completed. 
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